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Section I 


ORGANIZATION AND MANAGEMENT 

The 1985 Old Dominion University-NASA-Langley Research Center Summer Faculty 
Research Program, the twenty-second such institute to be held at the Langley 
Research Center, was planned by a committee consisting of the Co-Directors 
and Langley staff members from the Research Divisions and the Office of 
University Affairs. 

Each individual applying for the program was provided a listing of research 
problems available to the Langley Fellows. Each individual was requested to 
indicate his or her problem preference by letter to the University Co-Director. 
The desire to provide each Fellow with a research project to his/her liking 
was given serious consideration. 

An initial assessment of the applicant's credentials was made by the University 
Co-Director. The purpose of this assessment was to ascertain to which Divisions 
the applicant's credentials should oe circulated for review. Each application 
was then annotated reflecting the Division to which the applications should be 
circulated. On February 27 after the applications had been reviewed by the 
various divisions, a committee consisting of staff members from the various 
Divisions, the University Affairs Officer and the University Co-Director met. 

At this meeting the representatives from the various Divisions indicated those 
individuals selected by the Divisions. 

The University Co-Director then contacted each selected Fellow by phone extending 
the individual the appointment. The University Co-Director also forwarded each 
selected Fellow a formal letter of appointment confirming the phone call. 
Individuals were given ten days to respond in writing to the appointment- After 
letters of acceptance were in hand, the NASA-Langley University Affairs Officer 
contacted the various Directorate Technical Assistants advising them of their 
Fellows for the summer program. 

Each Fellow accepting the appointment was provided material relevant to housing, 
travel, payroll distribution and a listing of all NASA-Langley Research Fellows. 
Each Fellow, in advance of commencing the program, was contacted by his or her 
Research Associate or a representative of the branch. 

At the assembly meeting, the host, the NASA-Langley University Affairs Officer, 
introduced the Langley Research Center Chief Scientist, who formally welcomed 
the Summer Fellows. Mrs. Jane Hess, Head, Technical Library Branch, then 
briefed the Fellows on the use of the library. Mr. John Banks, Manager, Langley 
Research Cafeteria, briefed the Fellows relevant to the cafeteria policies, 
hours etc. Mr. Roger Butler of the Computer Management Branch briefed the 
Fellows on the Computational Facilities. The subject of security at the Langley 
Research Center was discussed by 0. J. Cole, Chief, Security Branch. 

Throughout the program the University Co-Director served as the principal liaison 
person and haa frequent contacts with the Fellows. The University Co-Director 
also served as the principal administrative officer. At the conclusion of the 
program, each Fellow submitted an abstract describing his/her accomplishments. 



Each Fellow was also required to conduct a fifteen minute oral presentation of 
his/her accomplishments as part of the lecture series. Each Fellow was re- 
quested to complete a questionnaire provided him/her for purposes of evaluation 
of the Summer Program. 
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Section II 


RECRUITMENT AND SELECTION OF FELLOWS 


Returning Fellows 

An invitation to apply and participate in the Old Dominion University-Langley 
Research Center Program was extended those individuals who held previous 
Langley Fellow appointments. Twenty-four individuals responded to the invi- 
tation; however, only eighteen were selected. Although eighteen previous Langley 
Fellows participated, seven individuals, who were participants in other similar 
programs, also participated in the program bringing the total number of returnees 
to twenty-five. 


New Fellows 


Although ASEE distributed a combined brochure of the Summer Programs, more than two 
hundred personal letters were mailed to Deans and Department Heads of various 
engineering schools in the East, South and Midwest requesting their assistance in 
bringing to the attention of their faculties the Old Dominion University-Langley 
Research Center Program. In addition to the above, a number of departments of 
chemistry, physics, computer science and mathematics at colleges (including 
community colleges) and universities in the State of Virginia as well as neighboring 
states were contacted regarding this program. Although minority schools in 
Virginia and neighboring states were included in the mailing, the alternating Co- 
Di rector from Hampton University made site visits to minority schools soliciting 
applicants. These efforts resulted in a total of one hundred eighteen formal 
applications, all indicating the Old Dominion University-Langley Research Center 
Program was their first choice and a total of fifty indicating ODU-LaRC Program 
was their second choice. The total number of applications received came to one 
hundred sixty-eight (Table 1). 

Forty-five applicants formally accepted the invitation to participate in the program. 
The average age of the participants was 41.3. On June 1, 1985, the University 
Co-Director received a letter from one of the appointees informing us that he was 
withdrawing from participating in the program. A replacement was available and 
the appointment made maintaining the total of forty-five participants in the 
program. 
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TABLE I 


First Choice Applications 


Total 

Females 

Males 

Minority Schools 


Black 

NonBlack 

I 

Black ! 

1 

NonBlack 

Represented 

118 

3 

4 

is i 

96 

12 



Second 

Choice Applications 


Total 

Females 

Males 

Minority Schools 


Black 

1 

! NonBlack 

Black 

NonBlack 

Represented 

50 

! 

0 

1 

5 

3 

42 

6 


NASA-LaRC FELLOWS 


TOTAL 

Females 

Mai 

es 

Minority Schools 


Black 

NonBlack 

Black 

NonBlack 

Represented 

45 

2 

2 

5 

36 

4 


First Year Fellows 

Returnees 



20 


25 
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Section III 


STIPENDS AND TRAVEL 

A ten-week stipend of $6,500 was awarded to each Fellow. Although this stipend 
was a considerable improvement over previous years, it still fell short (for 
the majority of Fellows) of matching what they could have earned based on their 
university academic salaries. This decision on their part does, however, 
clearly reflect the willingness of the Fellow to make some financial sacrifice 
in order to participate in the Summer Program. 

Travel expenses incurred by the Fellows from their homes to Hampton, Virginia, 
and return were reimbursed in accordance with current State regulations. 

Section IV 

LECTURE SERIES, TOURS, AND PICNIC 
Lecture Series 


Two lecture series were arranged for the NASA-ASEE Langley Fellows. The first 
series arranged consisted of general interest and specifically centered about 
research activities at the Langley Research Center. The second lecture series 
consisted of a number of lectures in Numerical Aerodynamics, specifically in 
Transition and Turbulence Modeling. These same lectures were part of a short 
course in Numerical Aerodynamics conducted by Old Dominion University's Mechan- 
ical Engineering Departments' Institute of Computational and Applied Mechanics. 
Although the Fellows were invited to all lectures, the following three lectures 
were for their benefit: Two lectures on Basic Statistical Mechanics of Turbu- 

lence and Statistical Theories of Turbulence and one lecture on Turbulent Shear 
Flows-Physics, Phenomenology and Control. The intent behind the Numerical 
Aerodynamics lecture series was to provide the Fellows with an exposure to 
Computational Fluid Dynamics and how valuable a tool it is in research at Langley. 

A complete listing of both lecture series is included in appendix III of this 
report. 


Tours and Picnic 


A briefing and tour of the NASA Langley Research Center was arranged for June 4, 
1985. A briefing and tour of the NASA Computational Facilities was arranged 
for June 5, 1985. A tour of the Acoustics Laboratory was arranged for July 23, 
1985, and a tour of the Full-Scale Wind Tunnel was arranged for July 26, 1985. 

A picnic for the Fellows, their families, and guests was held on June 14, 1985. 
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Section V 


RESEARCH PARTICIPATION 


The 1985 Ola Dominion Um vers ity-Lang ley Research Program, as in the past years, 
placed greatest emphasis on the research aspects of the program. Included in 
this report are abstracts from the Fellows showing their accomplishments during 
the summer. These abstracts, together with the comments of the Langley Research 
Associates with whom the Fellows worked, provide convincing evidence of the con- 
tinued success of this part of the program. The Fellows' comments during the 
evaluation of the program indicated their satisfaction with their research pro- 
jects as well as with the facilities available to them. 

The research projects undertaken by the Fellows were greatly diversified as is 
reflected in their Summer Research Assignments. Their assignments were as follows 


Number of Fellows 
Assigned 

2 

4 
6 
2 

5 
1 
1 
4 
2 
4 
1 

3 

4 
1 
1 
1 
1 
1 
1 


Division 


Analysis and Computational Division 
Instrument Research Division 
Flight Dynamics and Control Division 
Flight Electronics Division 
Flight Control Systems Division 
Materials Division 
Acoustics Division 
Structures and Dynamics Division 
Loads and Aeroelasticity Division 
Transonic Aerodynamics Division 
High-Speed Aerodynamics Division 
Atmospheric Sciences Division 
Space Systems Division 
Systems Engineering Division 
Facilities Engineering Division 
Personnel Division 
Management Support Division 
Financial Management Division 
Low-Speed Aerodynamics Division 


All but six of the participants were holders of the doctorate degree. The group 
was a highly diversified one with respect to background. In terms of their last 
degree, the breakdown revealed: 


Number 


Last Degree 


4 

1 

4 

1 

1 

1 

1 

1 

1 


Aeronautical Engineering 

Aerospace Engineering 

Chemistry 

Civil Engineering 

Computer Science 

Education 

Accounting 

Nuclear Chemistry 

Language and Literature 
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1 

1 

1 

1 

3 

2 

1 

6 

7 

1 

1 

3 

1 

1 


Mathematics/Computer Science 
Psychology 

Pharmaceutical Chemistry 
History 

Electrical Engineering 
Engineering Mechanics 
Astronomy 
Mathematics 

Mechanical Engineering 

Business Administration 

Systems Engineering/Computer Sciences 

Physics 

Operations Research 
Economics 


Extensions 


A portion of the funds that remained in the travel budget was used to grant 
a one week extension to fourteen Fellows in the program. For an extension to be 
granted approval by both the Research Associate and the University Co-Director 
was necessary. The following individuals were granted a one week extension: 

A. Blayton 

K. Brown 
R. Causey 
R. Coleman 
A. Dalessio 
M. Doria 

M. Ferguson 
R. Kapania 
J. Liburdy 
D. Lindner 
H. Robertshaw 
G. Rublein 

L. Silverberg 
P. Skoner 


Attendance of Short Courses, Seminars and Conferences 


During the course of the summer there were a number of short courses, seminars 
and conferences, the subject matter of which had relevance to Fellows' research 
projects. A number of Fellows requested approval to attend one or more of these 
conferences as it was their considered opinion that the knowledge gained by their 
attendance would be of value to their research projects. Those Fellows who did 
attend had the approval of both the Research Associate and the University Co- 
Di rector. The following is a listing of those Fellows attending either a short 
course, seminar or conference. 
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C. Cox attended a SIAM meeting in Pittsburg on June 24-26. He also presented 
a paper entitled, "On the Accuracy of Least-Squares Methods in the Presence of 
Corner Singularities." 

M. Maughmer attended a workshop held at the University of Notre Dame on 
June 17-19. 

G. Selby attended a short course on Fluid Mechanics Measurements held at the 
University of Minnesota on June 17-21. G. Selby's employer. Old Dominion Univer- 
sity, provided all funds necessary for the trip and course. 

F. Gross attended an IEEE/APS and URSI Symposium held in Vancouver on 
June 17-21. F. Gross represented the Antenna and Microwave Research Branch, FED 
at this symposium. His Research Associate requested that the ASEE program pro- 
vide reimbursement for certain of his expenses. F. Gross was reimbursed for 
$115 in accordance with his request. 

E. Lawrence attended Texas Instruments Industrial and Development Center in 
Johnston City, Tennessee, on July 8-19 and he also attended the Allen-Bradley 
Company's Training Center in Philadelphia, Pennsylvania, on July 22-27. In both 
cases the purpose of the visit was to allow the Fellow to attend short courses 
on Programmable Controllers (PC's) Programming. Both companies provided tuition- 
free enrollment for the Fellow. 

R. Longman and H. Robertshaw both attended the VPI and SU/AIAA Symposium 
on Control of Large Structures held in Blacksburg, Virginia, on June 12-14. 


Anticipated Publications Resulting from 
” Fellows Research Efforts 


Turbulence for Flight - to be submitted to Journal of Aircraft - 
G. Trevino. 

System Indentification Using a Recursive Form of the ERA 
Method - R. Longman. 

Vapor Screening in a Cryogenic Environment - to be submitted to 
AIAA Journal - G. Selby. 

Design of a Variable Coefficient Viscous Inertial Damper Using 
the Dual -General i zed Hessenberg Representation - to be 
submitted to Alberton Conference - D. Lindner. 

Phase of Gain Margins in Grassman Manifolds - to be submitted to 
IEEE Transactions - G. Rublein. 

Effects of Radiation on a Polyarylene Ether - R. Kiefer. 

Aspects of Resolution Concerning Multi-particle Dectection Using 
Far-Field Holography - to be submitted to Applied Optics - 
J. Liburdy. 
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A Control System Design Approach for Flexible Spacecraft - to be 
submitted to AIAA Journal of Spacecraft and Rockets - L. Silverberg. 

Putting Offline Documentation Online: A Staged Approach - to be 

submitted to AMC - SI6D0C - P. Carlson. 


An Annotated Bibliography of Online Documentation - to be submitted 
to AMC - SI6D0C - P. Carlson. 


Equivalent Continuum Analysis of Lattice Structures - J. Dow. 

OPTIMAL: A Parse Tree Manipulation Language - to be submitted to 

IEEE Transactions or ACM Transactions - L. Morel 1. 

Atmospheric Profiles in Atmospheric Environments - R. Coleman. 

In addition to the above, six additional Fellows will be publishing their 
results . 


Anticipated Research Proposals 

Flow Visualization Studies in a Cryogenic Environment - to be 
submitted to NASA-Langley - .G. Selby. 

An Automated, Open Air, Absorber Measurement System - to be 
submitted to NASA-Langley - F. Gross. 

Design Application of Factored Representations of Linear Systems - 
to be submitted to NASA-Langley - G. Rublein. 

Radiation Effects on High-Performance Polymers - to be submitted 
to NASA-Langley - R. Kiefer. 

Particle Detection Using Far-Field Holography in Turbulent Gas 
Flows - to be submitted to NASA-Langley - J. Liburdy. 

Characterization of the Interaction Between the Non-Linear Lift- 
body Motion and the Linear Elastic Motion of Rotating Flexible 
Spacecraft - to be submitted to NASA-Langley - L. Silverberg. 

Experimental Determination of the Stiffness Matrices for Symmetric 
Laminates - to be submitted to NASA-Langley - H. Patrick. 

Self-Generating Manuals: A Way of Incorporating User Cotranent in 

Online Documentation - to be submitted to NASA-Langley - 
P. Carlson, L. Morell. 

Sodium Vapor Solar Powered Lasant - to be submitted to NASA-Langley - 
C. Blount. 
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OPTIMAL: A Parse Tree Manipulation Language - to be submitted to 

NASA-Langley - L. Morell. 

Identification Techniques for MAST - to be submitted to NASA-Langley 

- J. Lilly. 

Influence of Antennas Flexibility in the Design of COFS Experiment 

- to be submitted to NASA-Langley - R. Kapania. 

In addition to the above, ten additional Fellows anticipated submitting pro- 
posals either to NASA, NSF or other agencies. 
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Section VI 


SUMMARY OF PROGRAM EVALUATION 


A program evaluation questionnaire was provided each Fellow and each Research 
Associate involved with the program. A sample of each questionnaire appears in 
Appendix V of the report. Forty- two of the forty-five Fellows responded. Twenty- 
four of the forty-five Research Associates responded. The results of the question- 
naires are summarized below. 


Fe 1 1 ows 


Ninety-five percent of the Fellows indicated they had found adequate information 
relevant to the program. 

Eighty-eight percent of the Fellows indicated that they had an early contact 
with the anticipated Research Associate. Without exception the Fellows considered 
the early contact as essential. 

Eighty-eight percent of the Fellows indicated that they were given a choice 
of research projects. 

Ninety-eight percent of the Fellows indicated that the research project was 
both challenging and in their field of interest. 

Ninety percent of the Fellows indicated that they anticipate continuing with 
their research upon return to their institutions. 

Ninety-one percent of the first year Fellows look forward to returning to 
Langley for a second year should they be extended that invitation. 

Ninety- two percent of the Fellows indicated satisfaction with the lecture 
series. 

All but six of the Fellows considered the stipend to be adequate. One 
recommended the stipend be increased by $1,000. 

Eighty percent of the Fellows indicated that they had no difficulty finding 
housing. Two individuals indicated dissatisfaction with the accommodations at 
Hampton University. 

All Fellows indicated satisfaction with the administration of the program. 

Fourteen Fellows provided the titles of papers they anticipate publishing 
shortly. Six others have plans to also publish a paper. 

Twelve Fellows indicated that they definitely will submit proposals to NASA 
Langley to continue their research. Ten other Fellows are comtemplating doing so. 
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Fellows' Commen ts 


The comments were as follows: excellent program for teachers who come from 

schools where there is little opportunity for research; very impressed with pro- 
gram; supremely valuable; very worthwhile and extremely well administered; very 
useful for exposing University professors to engineering problems of current 
interest; believe NASA benefits from program by the infusion of new people for a 
brief period; participants benefit by association with talented individuals; 
environment at NASA Langley very conducive to doing research; lovely program, it 
seems to benefit everyone; excellent program, excellent opportunity for a young 
professor; administrators served us well, very well administered, administration 
is especially to be commended ; everything provided on time and with no administra- 
tive problems; everyone was courteous, helpful and friendly; very good management, 
couldn't be better. 


Fellows' Recommendations 


Recommendations included the following: a housing officer is needed to assist 

in locating housing for the Fellows; recommend some monies such as $l,500-$3,000 
be set aside to aid Fellows in obtaining necessary supplies; monies should be pro- 
vided for those research projects requiring travel; computer system shortcomings: 
(1) inadequate diagnostics, (2) inexperienced consulting staff, (3) 24-48 hour 
job turnaround time, (4) long waits for response at terminal coupled with a timed 
shutdown that doesn't permit one to go away and come back, (5) file editing can 
be impossible, (6) limited times of access to computer - evening and weekend use 
would help, (7) REWIND exists on no other system - get rid of it. 


Research Associates 


All Research Associates responding indicated that their Fellows were adequately 
prepared to undertake his/her research assignment. 

All Research Associates responding were satisfied with the diligence, interest 
and enthusiasm that their Fellows exhibited. Research Associates responding were 
satisfied with the productivity and results of their Fellows' summer effort. 

All Research Associates responding expressed an interest in serving in the 
program again. 


Research Associates' Comments and Recommendations 


Comments and recommendations are given below. 

The interaction with highly capable people and a "new look" at some problems 
seems highly useful. This program is an excellent means to provide LaRC with 
access to skilled researchers with fresh viewpoints. The Fellows appear to 
benefit as well because it gives them access to current research problems. With 
the right pairing this program is beneficial to both parties. This program is a 
good experience for both the Fellow and his Associate, I personally have benefited 
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greatly from this association. This program is a good idea. An increase in the 
number of fellowships would help. Also, it would help if the restriction to 
U.S. citizens could be dropped. This is an excellent program. When Fellows 
have been selected, notified and have accepted, then the NASA associate should 
be specifically contacted by the University Co-Director's office and asked to 
communicate with the researcher. This year it seemed a bit awkward. I wasn't 
sure if the offer was accepted etc., and wished to communicate in advance. Some 
directions to the NASA associate would be appreciated. What type of problems 
do the fellows want, how much interaction, and what direction? 
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Section VII 


CONCLUSIONS AND RECOMMENDATIONS 


Conclusions 


Comments advanced by the Research Fellows, as well as others associated with 
the program, indicated satisfaction with the program. 

All Fellows considered contact with the anticipated research associate, prior 
to arrival, as essential and most desirable. 

A majority of the Fellows considered the stipend as adequate; however, there 
were a few who would only agree to this after considering the peripheral benefits 
associated with the program. 

The overall quality of the Fellows was considered outstanding to excellent 
by all those associated with the program. 

All Fellows indicated the program should be continued. 


Recommendations 


The Fellows recommended the practice of contact by the research associate, 
prior to arrival, be continued. They further recommended a site visit to Langley, 
in advance of the commencement of the program, be favorably considered. 

The Fellows recommended the stipend be increased, one suggested by as much as 
$1,000. They indicated that some very qualified individuals are discouraged 
from applying with the current stipend. 

The Fellows considered the tenure of the program as too short and recommended 
it be lengthened to perhaps twelve weeks. 

Some research associates recommended recruiting more candidates like the ones 
they had this year. One research associate recommended allowing anticipated 
returnees a short visit to Langley during winter academic break to plan second 
summers' work. One research associate indicated it desirable not to assign more 
than on Fellow to an associate. 

All individuals associated with the program recommended that the program be 
conti nued. 
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1985 

NASA-ASEE FELLOWS 
ODU-NASA LANGLEY 
RETURNEES 


Fellow 

Age 

Assigned to 

Research Associate 

Mrs. Blayton, Angela 
Assistant Professor 
Business Department 
Hampton University 
Hampton, VA 23668 

38 

Financial Management Division 

Joseph C. Struhar 
Building 1192 
Tel. 865-2724 

Dr. Blount, Charles E. 
Professor 

Physics Department 
Texas Christian University 
Fort Worth, TX 76129 

54 

Space Systems Division 

Nelson W. Jalufka 
Building 1247A 
Tel. 865-3781 

Dr. Brandt, Daniel A. 

Professor 

Industrial Technology Department 
Western Wisconsin Technical Institute 
Lacrosse, WI 54601 

53 

Flight Dynamics and Control Division 

Marion A. Wise 
Building 1268A 
Tel. 865-3871 

Dr. Brown, Kenneth G. 
Associate Professor 
Chemical Science Department 
Old Dominion University 
Norfolk, VA 23508 

41 

Instrument Research Division 

George M. Wood 
Building 1230 
Tel. 865-2466 

Dr. Bynum, William L. 
Professor 

49 

Flight Dynamics and Control Division 

Fenton W. Harrison 
Building 1268A 

Computer Science Department 
College of William and Mary 
Williamsburg, VA 23185 


- 

Tel. 865-3871 



40 


Dr. Carlson, Patricia A. 

Associate Professor 
Humanities Department 
Rose-Hulman Institute of Technology 
Terre Haute, IN 47803 

Dr. Causey, Robert L. 

Professor 

Computer Science Department 
Christopher Newport College 
Newport News, VA 23606 

Dr. Coleman, Randolph A. 

Associate Professor 
Chemistry Department 
College of William and Mary 
Williamsburg, VA 23185 

Dr. Dalessio, Anthony T. 

Assistant Professor 
Psychology Department 
Old Dominion University 
Norfolk, VA 23508 

Dr. Doria, Michael L. 

Associate Professor 
Mechanical Engineering Department 
Valparaiso University 
Valparaiso, IN 46383 

Dr. Dow, John 0. 

Assistant Professor 
Civil Engineering Department 
University of Colorado 
Boulder, CO 80309 


58 


41 


33 


46 


44 


Analysis and Computation Division 


Flight Control Systems Division 


Atmospheric Sciences Division 


Personnel Division 


Transonic Aerodynamics Division 


Structures and Dynamics Division 


Susan J. Voigt 
Building 1221 
Tel. 865-2083 


Ricky W. Butler 
Building 1220 
Tel. 865-3681 


W. Cofer/R. Harriss 
Building 1202 
Tel. 865-4372/3237 


L. Cornell Burcher 
Building 1219 
Tel. 865-2165 


Manuel D. Salas 
Building 1244 
Tel. 865-2627 


Harold G. Bush 
Building 1148 
Tel. 865-2498 



00 


Dr. Gosink, Thomas A. 
Research Associate 
Geochemistry Department 
Geophysical Institute 
University of Alaska 
Fairbanks, AK 99701 

49 

Atmospheric Sciences Division 

Leonard McMaster 
Building 1250 
Tel. 865-2065 

Dr. Kandil, Osama A. 
Associate Professor 
Mechanical Engineering and 
Mechanics Department 
Old Dominion University 
Norfolk, VA 23508 

41 

Transonic Aerodynamics Division 

Manuel D. Salas 
Building 1244 
Tel. 865-2627 

Dr. Kyrala, All 
Professor 
Physics Department 
Arizona State University 
Tempe, AZ 85282 

64 

Flight Control Systems Division 

Felix L. Pitts 
Building 1220 
Tel. 865-3681 

Dr. Liburdy, James A. 

Associate Professor 
Mechanical Engineering Department 
Clemson University 
Clemson, SC 29631 

36 

High-Speed Aerodynamics Division 

Dennis Bushnell 
Building 1247A 
Tel. 865-4546 

Dr. Lilly, John H. 

Assistant Professor 
Electrical Engineering Department 
University of Kentucky 
Lexington, KY 40506 

36 

Flight Dynamics and Control Division 

Suresh M. Joshi 
Building 1268A 
Tel. 865-4591 

Dr. Longman, Richard W. 

Professor 

Mechanical Engineering Department 
Columbia University 
New York, NY 10027 

42 

Structures and Dynamics Division 

Jer-Nan Juang 
Building 1293B 
Tel. 865-2881 



Dr. Mattix, Larry 
Associate Professor 
Chemistry/Physics Department 
Norfolk State University 
Norfolk, VA 23504 

38 

• 

Instrument Research Division 

William Winfree 
Building 1230 
Tel. 865-3036 

Dr. Nebel , Henry 
Associate Professor 
Physics Department 
Alfred University 
Alfred, NY 14802 

45 

Atmospheric Sciences Division 

James M. Russell 
Building 1250 
Tel. 865-4789 

Mr. Nold, Dean E. 

Professor 

Electrical Engineering Technology 
Department 

Purdue University Calumet 
Hammond, IN 46323 

53 

Flight Control Systems Division 

Jack J. Hatfield 
Building 1298 
Tel. 865-2171 

Dr. Robertshaw, Harry H. 

Associate Professor 
Mechanical Engineering Department 
Virginia Polytechnic Institute 
and State University 
Blacksburg, VA 24061 

43 

Structures and Dynamics Division 

Garnett C. Horner 
Building 1293B 
Tel. 865-2817 

Dr. Sanwal, Jagdish C. 
Associate Professor 
Mathematics Department 
College of William and Mary 
Williamsburg, VA 23185 

51 

Flight Dynamics and Control Division 

Nancy E. Orlando 
Building 1268A 
Tel. 865-3871 

Dr. Scanlon, Charles H. 
Associate Professor 
Computer Science, Math and 
Physics Department 
Arkansas State University 
State University, AR 72401 

48 

Flight Control Systems Division 

Charles E. Knox 
Building 1168 
Tel. 865-3621 



ro 

o 


Dr. Selby, Gregory V. 

36 

Transonic Aerodynamics Division 

Richard Barnwell 

Assistant Professor 


Building 1229 

Mechanical Engineering and 
Mechanics Department 
Old Dominion University 
Norfolk, VA 23508 



Tel. 865-2601 

Dr. Trevino, George 

43 

Flight Control Systems Division 

Roland Bowles 

Associate Professor 


Building 1168 

Mechanical Engineering and 
Engineering Mechanics Department 
Michigan Technological University 
Houghton, MI 49931 



Tel. 865-3621 
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1985 

NASA-ASEE FELLOWS 
ODU-NASA LANGLEY 
1st YEAR 


Fellow 

Age 

Assigned To 

Research Associate 

Dr. Chachere, Bernadette P. 
Chairman 

Economics Department 
Hampton University 
Hampton, VA 23668 

41 

Flight Electronics Division 

Glenn R. Taylor 
Building 1299 
Tel. 865-3541 

Dr. Chen, Kwan-Yu 
Professor 

Astronomy Department 
University of Florida 
Gainesville, FL 32611 

55 

Space Systems Division 

Edmund J. Conway 
Building 1247A 
Tel. 865-3781 

Dr. Cox, Christopher L. 
Assistant Professor 
Mathematics Department 
Clemson University 
Clemson, SC 29631 

29 

Loads and Aeroelasticity Division 

John W. Edwards 
Building 648 
Tel. 865-4236 

Dr. Day, Henry P. 

Assistant Professor 
Virginia Polytechnic Institute 
and State University 
Blacksburg, VA 240bl 

37 

Low-Speed Aerodynamics Division 

Bruce J. Holmes 
Building 1212 
Tel. 865-2877 

Mr. Ferguson, Milton W. 
Assistant Professor 
Physics Department 
Norfolk State University 
Norfolk, VA 23504 

47 

Instrument Research Division 

Sheila T. Long 
Building 1267 
Tel. 865-3892 



34 


Flight Electronics Division 


Dr. Gross, Frank B. 

Assistant Professor 
Electrical Engineering Department 
Florida State University 
Tallahassee, FL 32308 

Dr. Hansen, James R. 

Assistant Professor 
History Department 
University of Maine 
Orono, ME 04473 

Dr. Kapania, Rakesh K. 

Assistant Professor 
Aerospace and Ocean Engineering 
Department 

Virginia Polytechnic Institute 
and State University 
Blacksburg, VA 24061 

Dr. Kiefer, Richard L. 

Professor and Chairman 
Chemistry Department 
College of William and Mary 
Williamsburg, VA 23185 

Dr. Lawrence, Ellis E. 

Associate Professor 
Industrial Technology Department 
Elizabeth City State University 
Elizabeth City, NC 27909 

Dr. Leemis, Lawrence M. 

Assistant Professor 
Industrial Engineering Department 
University of Oklahoma 
Norman, OK 73069 


33 


Management Support Division 


29 Structures and Dynamics Division 


48 Materials Division 


36 Facilities Engineering Division 


28 Space Systems Division 


Marion C. Bailey 
Building 1299 
Tel. 865-3631 


Richard T. Layman 
Building 1151 
Tel. 865-3511 


Brantley R. Hanks 
Building 1293B 
Tel. 865-3055 


George F. Sykes 
Building 1205 
Tel. 865-4555 


John L. Gilbert 
Building 1209 
Tel. 865-4534 


W. Douglas Morris 
Building 1232 
Tel. 865-2768 


ro 

u> 



ro 
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Dr. Lindner, Douglas K. 

Assistant Professor 
Electrical Engineering Department 
Virginia Polytechnic Institute 
and State University 
Blacksburg, VA 24061 

33 

Flight Dynamics and Control Division 

R. C. Montgomery 
Building 1268A 
Tel. 865-4591 

Dr. Maughmer, Mark D. 

Assistant Professor 
Aerospace Engineering Department 
The Pennsylvania State University 
University Park, PA 16802 

35 

Transonic Aerodynamics Division 

William D. Harvey 
Building 641 
Tel. 865-2631 

Dr. Morel 1, Larry J. 
Assistant Professor 
Computer Science Department 
College of William and Mary 
Williamsburg, VA 23185 

33 

Analysis and Computation Division 

Edmond H. Senn 
Building 1268A 
Tel. 865-2558 

Dr. Pacheco, Maria D. 

Visiting Assistant Professor 
Science and Math Division 
College of The Virgin Islands 
St. Thomas, U.S.V.I. 00802 

28 

Instrument Research Division 

Jag J. Singh 
Building 1230 
Tel. 865-3907 

Mr. Patrick, Howard V. 

Assistant Professor 
Aeronautical Engineering Department 
Embry-Riddle Aeronautical University 
Regional Airport 
Daytona Beach, FL 32014 

49 

Acoustics Division 

John S. Mixson 
Building 1208 
Tel. 865-3561 

Dr. Rublein, George T. 
Associate Professor 
Mathematics Department 
College of William and Mary 
Williamsburg, VA 23185 

50 

Flight Dynamics and Control Division 

Douglas B. Price 
Building 1268A 
Tel. 865-4681 



J 


Dr. Silverberg, Larry M. 

Assistant Professor 
Mechanical & Aerospace Engineering 
Department 

North Carolina State University 

28 

Space Systems Division 

L. Bernard Garrett 
Building 1232 
Tel. 865-3666 

Mr. Skoner, Peter R. 
Instructor of Physics and 
Management 
St. Francis College 
Loretto, PA 15940 

28 

Systems Engineering Division 

Charles E. Cockrell 
Building 1209 
Tel. 865-4666 

Mr. Turner, James C. 
Graduate Teaching Assistant 
Mathematics Department 
Carnegie-Mel Ion University 
Pittsburg, PA 15213 

38 

Loads and Aeroelasticity Division 

Allan R. Wieting 
Building 1265 
Tel. 865-3423 


ro 

tn 
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ASEE-NASA 

Old Dominion University - Langley Research Center 
Lecture Series 
General 

Location: Building 1219, Room 225 

Time: 9:00 a.m. - 10:30 a.m. 


Date 

Speaker 

Topic 

June 7 

Mr. Robert E. Bower 
Associate Director 
NASA/Langley Research Center 

An Overview of the Langley 
Research Center 

June 19 

Dr. Douglas L. Dwoyer, Branch Head 
Computational Methods Branch 
High-Speed Aeronautics Division 

Computational Fluid Dynamics 

June 28 

Dr. Roger A. Breckenridge 
LaRC Space Station Office 

NASA's Space Station 

July 2 

Mr. Frederick 0. Allamby 
Space Flight Experiment Definition 
and Integration Office 

Commercialization of Space 

July 11 

Dr. James H. Starnes, Branch Head 
Structural Mechanics Branch 
Structures and Dynamics Division 

Composites in Aeronautics and 
Space 

July 25 

Dr. William P. Winfree 
Instrument Research Division 

Applications of Quantitative 
Physical Measurements 

August 1 

Dr. Samuel E. Massenberg 
University Affairs Office 
NASA/LaRC 

Preparation and Submission 
of Unsolicited Proposals 

August 7 
8:15 a.m. 

ASEE-NASA Fellows assigned to Aeronautics Directorate 

Dr. Henry Day Dr. Mark Maughmer 

Dr. Michael Doria Dr. James Liburdy 


ASEE-NASA Fellows assigned to Electronics Directorate 


Dr. 

Patricia Carlson 

Dr. 

Dr. 

Larry Morel 1 

Dr. 

Dr. 

Ali Kyra la 

Dr. 

Dr. 

Robert Causey 

Dr. 

Mr. 

Dean Nold 

Dr. 

Dr. 

Charles Scanlon 

Mr. 

Dr. 

John Lilly 

Dr. 

Dr. 

Dr. 

Douglas Lindner 
William Bynum 

Dr. 


Jagdish Sanwal 
George Rublein 
Frank Gross 
Bernadette Chachere 
Kenneth Brown 
Milton Ferguson 
Larry Mattix 
Maria Pacheco 
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August 8 
8: 15 a.m. 


ASEE-NASA Fellows assigned to Space Directorate 


Dr. Randolph Coleman Dr. Charles Blount 

Dr. Henry Nebel Dr. Kwan-Yu Chen 

Dr. Lawrence Leemis Dr. Thomas Gosink 

Dr. Larry Silverberg 

ASEE-NASA Fellows assigned to Structures Directorate 

Mr. Howard Patrick Dr. Harry Robertshaw 

Dr. Christopher Cox Dr. John Dow 

Mr. James Turner Dr. Richard Longman 

Dr. Richard Kiefer Dr. Rakesh Kapama 

ASEE-NASA Fellows assigned to Management and Operations 
Directorate 

Mrs. Angela Blayton Dr. James Hansen 

Dr. Anthony Dalessio 

ASEE-NASA Fellow assigned to Systems Engineering and 
Operations Directorate 

Dr. Ellis Lawrence 
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ASEE-NASA 

Old Dominion University - Langley Research Center 
Lecture Series 
Numerical Aerodynamics 

Location: Building 1244, Room 223 

Time: 10:00 A.M. - 12 Noon 


Date 

Speaker 

Topic 

June 24 

Dr. T. A. Zang 

NASA Langley Research Center 

Introduction to Spectral Methods 

June 25-27 

M. G. Worster 

Massachusetts Institute of 
Technology 

Introduction to Hydrodynamic 
Stability, Inviscid Instability 
of a Vortex Sheet (Free Shear 
Layer), and Two-Dimensional 
Parallel Shear Flow (Inviscid 
Disturbances) 

July 1 

M. G. Worster 

Massachusetts Institute of 
Technology 

Taylor-Couette Flow (Linear 
Theory) 

July 2-3 

Annick Pouquet 
Observatoire de Nice, France 

Basic Statistical Mechanics of 
Turbulence and Statistical 
Theories of Turbulence 

July 8 

D. M. Bushnell 

NASA Langley Research Center 

Turbulent Shear Flows -Physics , 
Phenomenology, and Control 

July 9-11 

M. G. Worster 

Massachusetts Institute of 
Technology 

Two-Dimensional Parallel Shear 
Flows (Viscous Disturbances), 
Linear and Nonlinear Theories 

July 15-18 

Tom Eidson 

Georgia Institute of Technology 

Inadequacies of Reynolds Averaging 
of Turbulence, Simulation of 
the Flow Equations and Subgrid 
Modeling Theory, and Advanced 
Topics and Problems in Large 
Eddy Simulation (LES) 

July 22 

Dr. T. A. Zang 

NASA Langley Research Center 

Numerical Methods for Simulation 
of Transition and Turbulence 

July 23-25 

Tom Eidson 

Georgia Institute of Technology 

Data Analysis of Large Eddy 
Simulation (LES), LES of 
Turbulent Convection, and LES 
of Channel Flow 
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A STUDY OF THE FAB RATE 


Angela Macklin Blayton 
Assistant Professor 
Accounting and Finance 
Hampton University 
Hampton, Virginia 


The Fabrication Division (Fab) is responsible for providing 
and procuring the services required to manufacture and test 
aeronautical and aerospace research hardware, related ground 
support equipment, and research facilities test equipment. In 
supporting the engineering and research organizations at the 
Center, Fab also performs experimental and developmental work 
related to the use of new materials processing techniques. 

Fab by itself cannot support all of the Center's fabrication 
requirements; therefore, one-half of the work is performed 
out-of-house by contractors. As a result, there was a lot of 
competition for Fab's resources which required that the 
researcher only pay for material costs. If the researcher had 
the same work performed out-of-house, he would have to pay the 
full cost of materials and labor; therefore, in light of the 
competition for in-house work, the Center developed a "Standard 
Rate" which made the hourly cost of the job the same irrespective 
to whether or not the job was performed in-house or out-of-house. 

The Fab Rate System is an operational cost recovery system. 
All stock issued, PR/PO purchased supplies, materials, equipment 
replacement and repair costs plus small business tasks and other 
contract costs must be recovered via the Fab System. The Fab 
rate is determined by dividing Total Operational Cost by Hours of 
Research Support and the Fabrication Division has the objective 
of delivering the hours needed to match recorded costs in an 
effort to maintain the predetermined Fab rate. 

The Program and Resources Division (PRD) is concerned with 
monitoring the Fab rate to insure its reasonableness through 
year-end. PRD is also concerned with controlling the dollar 
amount of uncosted obligations at year-end. 

The Financial Management Division ( FMD) is concerned with 
maintaining accurate records of actual costs along with assisting 
Fab and PRD with controlling the Fab rate, and as a represent- 
ative of FMD, I have met with Jeannie Duncan (PRD) and Bob Magee 
(Fab) at various intervals throughout the summer to discuss the 
status of the Fab rate which was set at $15.25 for Fiscal Year, 
1985. 



At our last meeting, I presented an analysis of Fab 
Contracts and a year-end cost projection for FY '85. Bob Magee 
presented a year-end labor hour projection for FY *85 along with 
an estimate of new Fab dollars to be obligated before 9/30/85. 


1. $9,672,743 

2. $8,555,191 

3. $649 , OOOhr 

4. $1,086,000 

5. $ 636,000 

6. $ 412,500 

7. $ 400,000 


FMD Year-End Costs Projection (see Exhibit 
A) 

Fiscal Year To Date Obligations as of 
6/30/85 

Fab's Year-End Projection 

Fab's estimated increase in Task Contract 
Obligations prior to Year-End 1985 

The additional dollars which should be 
costed ($1,086,000-450,000). 

Fab's estimated increase in obligations for 
support items. 

Obligation for Equipment (not on the books) 
Projections 


Obligations 
2. $8,555,191 

4. 1,086,000 

6. 412,500 

7. 400,000 
10,453,691 

1,509,911 FYTD Carryover 
11,963,602 Total Obligations 


Costs 

1. 9,672,743 

5. 636,000 

10,308,743 total 
costs 


Hours 

649,000 

Uncosted 

Obligations Cost Obi igations 

$11,963,602 $10,308,743 $1,654,859 

Projected Rate at Year-End - $ 10 , 308 ,743 = $15.89 

649,000 

Projected Rate - $15.89 

Predetermined Current Rate - 15.25 

Difference - .64 
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Recommendations : 


. Do not change the rate at this time. 

. Attempt to reduce the Fab hours put into 
AJO's and increase the RJO effort, cut back 
on stock issues, and put a lid on new obli- 
gations. 

. Continue discussions and analyses. 



Exhibit A 

An Analysis of FAB Contracts 
For The Month Ending June 30, 1985 


Contract # 
NAS-1-15416 

NAS-1-15876 

NAS-1-16369 

NAS-1-16610 


NAS-1-16613 

NAS-1-17222 

NAS-1-17514 

NAS-1-17986 


Other Direct Cost and Obligations 


Not Included Above: 
M1305 

$ 676,680 

$ 676,680 

-0- 

M1306 

1,487,385 

1,487,385 

-0- 

Ml 585 

1,392,377 

1,392,377 

-0- 

M1713 

58 

58 

-0- 

M1800 (Adj • ) 

(1,807,237) 

(1,807,237) 

-0- 

Subtotal 

1,749,263 

1,749,263 

-0- 

Total Direct 

$26,065,382 

$24,221,45$ 

$1,784,426 


Job Order t 
Ml 7 1 3 , Ml 305 , 
Ml 777 

Ml 305 , Ml 585 , 
M1773 

Ml 305 , M1711, 
M1775 

Ml 305 , Ml 71 4 , 
M1752 , Ml 765 , 
Ml 778 

M1712, Ml 776 
Ml 306 , M1766 
M1820 
M1773 

Subtotal 


PYTD 

Obi igat ions 
$6,696,169 

1,661,595 

5,456,474 

1,442,562 


1,144,619 

3,355,200 

4,250,000 

250.000 


PYTD 

Accruals 

$6,561,726 

1,661,595 

5,456,454 

1,393,376 


950,253 

3,346,121 

2,940,794 

161,874 

$22,472.1?? 


PYTD 

Uncosted Obligations 
$134,443 


49,186 


194,366 
9,079 
1,309,206 
88,126 
^ 1 778 474-25 


Exhibit A (Concluded) 


1 YTD 


Indirect Costs. 

Obligations 

Model and Fabrication 
support (Small Business)- 
M1765, Ml 30 5 

$3,633,072 

Fabrication Supplies- 
M1753, Ml 304 

2,788,457 

Equipment Maintenance, 
Repair & Replacement- 
M1752, Ml 303 , M1779 

2,136,229 

Supplies Purchase Request- 
M1774, Ml 586 

2,147,645 

Balance Call Agreement- 
M1769, Ml 581 

240,917 

RCA Contract Supplies- 
M1768, Ml 4 9 3 

238,152 

Machine Tool 

Rehabil l tat ion-Ml772 , M1584 

12,019 

Naval Air Rework Facility- 
M1770, Ml 582 

759 

Total indirect 

11,197,250 

Total 

$37,203,132 


PYTD PYTD 

Accruals Uncosted Obligations 

$3,459,159 $173,913 

2,788,457 -0- 


1,780,891 

2,011,948 

193,117 

235,902 

12,019 

759 

10,482,25 2 
$34,703,708 


355,338 

135,697 

47,800 

2.250 

- 0 - 

- 0 - 

714,998 

$2,499,424 


FYTD Cost Accruals as of 6/30/85 
Potential Accruals at Year-End, 9/30/85: 

Uncosted Obligations $2,499,424 

Less: Costs not to be Accrued inFY 1985 

Dynamic Engi neermg/NAS-1-17514 ( 350,468) 

Potential Year End Costs 


$7,523,787 

2,148,956 

$ 9,672,743 
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THE ULTRAVIOLET AND VISIBLE EMISSION OF SODIUM VAPOR 
PRODUCED BY LASER AND XENON CONTINUUM 


Charles E. Blount 
Professor 

Department of Physics 
Texas Christian University 
Fort Worth, Texas 76129 

It has been proposed that the space station will have a solar pumped 
laser. Presently, all proposed solar pumped laser systems require 
additional waste heat management equipment. The development of a solar 
powered metal vapor laser would eliminate the need for the heat management 
equipment. The object of this research has been to obtain the fluorescence 
spectra of sodium vapor excited by light similar to that of the sun, 
seeking those bands in the fluorescence spectra that might produce laser 
emission. 

Atomic alkali molecules have long been an object of considerable 
spectroscopic interest. Sodium molecules due to their well known spectrum, 
high absorption and emission cross section, and ease of production are of 
particular interest for laser systems. Presently, in the literature a 
large number of laser bands in the ultrviolet (1), visible (2) and infrared 
(3, 4) regions of the spectrum have been reported for sodium vapor. 

However, in all of these a laser was used for excitation. 

Of the fluorescence bands of sodium vapor observed in the ultraviolet and 
visible region of the spectrum and their dependence on the vapor pressure, 
future research should be directed on the band at 436 nm for which gain 
measurements have been reported (5) and at 451 nm. 


1) C. Y. R. Wu and J. K. Chen, Opt. Commun. JjO, (1984) 317-319. 

2) C. N. Man-Pichot and Alain Bril let, IEEE J. Quantum Electron, QE-16, 
(1980) 1103-1108. 

3) B. Wei legehausen, W. Luhs, A. Topouzkhaman, and J. d' Incan, Appl. 
Phys. Lett. 43 (1983) 912-914. 

4) W. Muller and I. V. Hertel, Opt. Commun. 45 ^ (1983) 400-402. 

5) J. T. Bahns and W. C. Stwalley, Appl. Phys. Lett. ^4, (1984) 826-827. 
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MECHANICAL FEATURES OF SPACE ROBOTICS 


Daniel A. Brandt 
Professor 

Industrial Technology Division 
Western Wisconsin Technical Institute 
LaCrosse, Wisconsin 

and 


Mechanical Engineering Department 
Milwaukee School of Engineering 
Milwaukee, Wisconsin 


This project became a series of projects, each involving some aspect of the 
mechanical design and analysis, of robotic mechanisms and features which 
will be utilized in space. 

Following is a brief summary of four of the these design analyses: 

1. Design of Space Tools -For reasons of safety and convenience, it is 
desirable that tools that are utilized in the outer space environment be 
actuated by robots instead of by astronauts. A screwdriver (see 
Appendix A), pliers (see Appendix B), and forceps (see Appendix C) , 
were designed, drafted and prototypes were manufactured and assembled. 
Each is powered and gripped by the end -effector qripper of the PUMA 600 
robot at Langley Research Center. 

2. End-Effector Force Analysis -A force analysis was made of the forces in 
• the end-effector mechanism of the robot at Langley. 

This analysis related (a) the output torque of the end-effector motor, 
(b) the tangential output force of the worm, (c) the tangential force 
transmitted to the worm gear, and (d) the resulting gripper force at the 
actuation end of the mechanism (see Appendix D). These relations were 
found to be: 
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= the force of the gripper 
= the tangential gear force 
= the tangential worm force 
= the actuating torque of the motor 


The maximum gripper force was found to be 80.2 pounds, based on a 
maximum input torque of 60 ounce-inches. 


3. Velocity Analysis of Epicyclic Gear Train in Robot Arm. -Joint #2 of the 
robotic mechanism employed an epicyclic gear train (see Appendix E). 

The velocity ratios of robot arm "A" to intermediate gears "C-D" and 
to the input motor shaft were found to be: 

_ ioyr _ A/gftfe _ /o.^f 

Nc A/ e ft c0 ^(a/c+A/s) 

where n is the angular velocity of each respective component and "N" is 
the number of teeth of each respective gear. 


4. Inertia Analysis of Robot Arm Mechanisms . -It was necessary to translate 
the inertia of two of the robot arms to an equivalent moment of inertia 
at the input motor shaft, (see Appendices F4G). Results included the 
fol lowi ng: 
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APPENDIX A 




Note Most of this is similar to pliers design 

(especially the right side of the mechanism) 
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THE DESIGN OF A MASS SPECTROMETER TOR ATMOSPHERIC 
ANALYSIS ON A TETHERED SATELLITE AND ENTRY VEHICLE 

Ke nne t h G B r own 

Associate Professor 
Department of Chemical Sciences 
Old Dominion University 
Norfolk, Virginia 23-08 


The use of mass spectrometry to sample the near 
surface environment of an entry test vehicle or low altitude 
tethered satellite is evaluated For successful operation in 
the higher pressure regions that exist at lower altitudes the 
mass spectrometer system must contain an inlet that i ) is 
capable of providing a sufficient pressure drop to protect 
the ion source, 2 ) can admit sample to the mass spectrometer 

that is representative of the external environment, and 3 ) 

the sample removal will not significantly perturb the flow at 
the surface of the vehicle In addition the remainder of the 
mass spectrometer must be located some distance from the 
inlet in order to be within the protective shielding As a 
result it is necessary to consider the design of the magnet 
in order to achieve a magnetic field that will allow the 
source and detection equipment to be situated apart f om the 
magnet and at the same time ensure that the total system is 
1 i gh twe i gh t 

The proposed inlet system consists of an array of 
capillary openings (approximately 10^ / cro^ ) The openings 

in the individual plates are uniform and range in diameter 
from 5 4 for a plate that is coated to 9 6ym for an 
uncoated plate The plates under study are either 4 0 0 yu> o r 
i 0 0 0 ym in thickness A plate is mounted in a holder with the 
surface of the plate flush with the surface of the holder 
The holder with the plate in place is mounted at the inlet of 
a quadrupole mass spectrometer The entire system is capable 
of being evacuated to 10"^ t or r or lower The test gas is 
admitted to / the system to obtain a flow which is parallel to 
the surface of the plate The pressure upstream and 
downstream of the plate is then measured as well as the 
pressure inside the mass spectrometer anc the intensities of 
the individual peaks that are present in the particular 
s amp I e 

The test gas used m the initial evaluation of the 
inlet was a miiture of 7* CO 2 in air The upstream pressure 
ranged from 1 to 1000 micrometer of Hg which corresponds to a 
mean free path range at room temperature of from 1 to 1 cm 
at 298K The peak intensities due to N 2 , O 2 » and CO 2 were 
determined as a function of applied pressure To evaluate the 
flow the peak intensities for each molecule were plotted as 
a function of the ratio of the transparency, z » of the plate 
to y/\L where X is the mean free path and L is the length of 
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the capillary (1) A sample plot is shown in figure i with a 
break in the curve indicating the onset of molecular flow 
occurring at approximately 7 for all three species Future 
experiments will involve plates of different thicknesses, 
higher gas temperatures and gases with a greater disparity in 
molecular weight 

The magnet in a mass spectrometer is usually designed 
to produce a homogeneous magnetic field As a result the 
object and image focus occur at the magnet face limiting the 
flexibility of the instrument An inhomogeneous magnetic 
field will allow tne object and linage points to be located 
away from the face of the magnet as shown in figure 2 The 
image plane, for different masses, will now be at an angle to 
the edge of the magnet instead of the parallel plane which 
exists for a homogeneous field We are currently 
investigating the effect of downsizing the magnet and 
changing the sector angle from the 100 degree case shown in 
figure 2 <2) to angles of 90 degrees or less Ray tracing, 
such as shown in figure 2, will enable the determination of 
the overall design for the system 

1 Johnson, J C , Stair, A T and Pritchard, J L , J Appl 

Phys 3_7 1 55 1 - 1 550 ( 1 9 6 6 ) 

2 Whitehead, T W and White, F A , Nuclear Instruments and 

Methods 103 437-445 (1972) 
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Figure 2 


. Calculated ion trajectories for an ion of mass 16 in a 24 inch 
magnet with an inhomogeneity of .5. Acceleration voltage of 2000keV 
and a of 500 gauss. 




SIMULATION OF A PARALLEL JAW END EFFECTOR 


Bill Bynum 
Professor 

Department of Computer Science 
College of William and Mary 
Williamsburg, V A 23185 

One of the current research projects in the Automation and Technology Branch on 
remote teleoperator control of robotic manipulators is the Distributed Artificially 
Intelligent System for Interacting with the Environment (DAISIE) (cf (1,2,31), a system 
developed by Nancy Orlando that combines artificially intelligent, goal-structured robot 
planning algorithms with traditional control methods for mechanical equipment In the 
DAISIE system, communication between the strategic planner and the various sensors 
and devices of the system is of paramount importance. When a new device is 
incorporated into the system, the set of commands that the device will accept and the 
responses that it will return must be specified; then this command protocol has to be 
encoded in a form usable by the processor that controls the device. In the case of an end 
effector, the specification of the command protocol is straightforward, but conveying the 
protocol to the processor controlling the end effector consists of burning the 
communciation protocol program into the ROM of the processor Unfortunately, it is 
impossible to forsee every eventuality when specifying a command protocol, which leads 
to iterating the previously described two-step process several times in developing a 
satisfactory command structure. The associated eipense and effort of this process justify 
the use of software simulation to develop a command protocol for the end effector. 

My project for the fellowship period has been to develop a collection of programs 
that will simulate the parallel jaw end effector used in the Intelligent Systems Research 
Lab (Wise, M.A., NASA Langley], The design goal for this collection of programs is to 
provide a graphical display of the end effector that moves in approximately real time and 
executes concurrently with the compiling DAISIE program. The program is to execute on 
the VAX 1 1/750 in the Intelligent Systems Research Lab. 

Graphics 

The first step of the project was to develop a graphical representation of the end 
effector that moves in approximately real time A program was developed in F0RTRAN77 
using the DI3000 graphics package driving either a DEC VT240 graphics terminal or a DEC 
VT100 terminal fitted with a Retrograde Graphics board. Unfortunately, this hardware 
could not change the graphics display quickly enough to approximate real time movement 
of the end effector. 

The program was rewritten for the VS1 1 graphics package driving the DEC VS1 1 
terminal. The VS1 1 terminal has four separate memories, called "memory planes' . in 
which images can be stored. This hardware obtains a higher display rate by overlapping 
the storage of the image in one of the planes while another of the planes is being 
displayed. For the sake of simplicity of program structure, only two of the memory 
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planes are used in the graphical display of the end effector Movement of the end 
effector from one gap opening to another is portrayed by subdividing the total gap into 
smaller parts and shoving the end effector in each of the intermediate positions. The 
speed of the movement is directly related to the number of subdivisions chosen and can 
be varied from barely perceptible to almost instantaneous. 

Concurrency 

Concurrency is simulated through the use of the VAX/VMS system subroutines, 
accessible from F0RTRAN77. These routines provide concurrently executing processes the 
capabilty of communication through "mailboxes" and of synchronization through "event 
flags " Additionally, there are subroutines with which a process can "create," or initiate 
execution of a separate process 

Program Structure 

The strategic planner for the DAISIE system is written in the version of LISP on the 
ISRL VAX, Franz LISP. This LISP provides the facility to incorporate dynamically into the 
LISP environment separately compiled FORTRAN subroutines and functions, which allows 
the graphics and concurrency capabilites described above to be controlled by a LISP 
program. 

An imported FORTRAN subroutine is called to start the simulation. This subroutine 
initializes the graphics display and the interprocess mailboxes and event flag clusters and 
initiates the subroutine that generates the graphics. The initialization subroutine 
eiecutes to completion, but cannot terminate until the process that it created, the 
graphics process, terminates Control returns to the calling LISP program, allowing 
additional LISP functions to be invoked. 

Each of the commands to the end effector is encoded as a LISP function . Each of 
these functions transforms its command into the format that the graphics process expects 
and invokes an imported FORTRAN subroutine to place the command message into the 
command mailboi and set an associated event flag. The graphics process is watching this 
event flag, waiting for it to be set. When the graphics process notices that a command has 
been given, it retrieves the command from the mailbox. The move or mitialize-jaw 
commands are enqueued for subsequent action, whereas the quit, status request, and 
debug commands receive an immediate reply Move commands are enqueued because, in 
general, they require a longer time to complete than a status request. It is possible for 
the LISP driver program to transmit commands, particularly status requests, and receive 
replies during the execution of a move command 

In the case of the request for status information, a reply consists of placing an 
appropriately formatted command m the reply mailbox and setting the associated event 
flag, the response to the move command is similar. The response to the debug command 
is to enable or disable the output to the user s terminal of various debugging information. 
The quit command causes the release of the system mailboxes and event flag cluster, as 
well as the termination of the LISP driver function, the graphics program, and the 
initialization process that created it. 

Sketches of the end effector for several jaw openings are shown in Figure 1. 
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ELECTRONIC SOFTWARE SUPPORT: 

ONLINE DOCUMENTATION, TUTORIALS, HELP AND ERROR MESSAGES 


Patricia Ann Carlson 
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Introduction 


Electronic software support has matured from the early days where 
terse and often cryptic error messages constituted the only form of 
online assistance. Currently, specific categories of electronic user 
support include: command help and error assistance; prompting; online 

tutorials, and online documentation. Driven by the need to make 
systems usable by a diversified community and the awareness among 
hardware and software vendors that "user-friendliness” is a 
marketable concept, man/machine interfaces and online assistance are 
central issues in research and development today. This report deals 
specifically with online and offline documentation. 

NASA/Lanqley Project 

I was asked by ay NASA colleagues — Dr. John Shoosmith and Susan J. 
Voigt, of the Analysis and Computation Division — to examine three 
issues — with particular attention to the Control Data Corporation 
environment and the NOS operating system: 

o the Langley users' documentation needs 

o the current library of paper documentation 

o the current capabilities of electronic computer 
support -- both at LaRC and state of the art 


Specific recommendations for implementing online user support at LaRC 
follow from the investigation of these three items. 

Overview of Findings 

Except for an influx of summer people, LaRC has a large but 
relatively fixed user community. The majority of NOS users are 
computer professionals or scientists/engineers working on specific 
analytical problems. Because the users are "computer literate," 
training and support typically begin at a relatively high level. For 
the most part users refer to the documentation for task-related 
Information — either procedural instructions to perform a function or 
explanat ions /de script ions of particular features. This report 
elaborates on the need for (1) easy-to-use, (2) timely, (3) accurate, 
(4) accessible, and (5) readable documentation in the LaRC user 
community. 
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Paper documentation is the traditional method for handling user 
information needs. Even so, only recently has serious attention been 
given to techniques for producing good offline documentation. Within 
the past three years, LaRC has identified a set of manuals that fit 
together in a logical and sequential fashion without unnecessary 
overlap and has set up comprehensive policies and procedures for 
writing and distributing quality paper documentation. This report 
addresses both the strengths (e.g., rigorous editorial and review 
processes) and the weaknesses (e.g., insufficient field testing and 
lack of modern document design techniques) of this system. 

The efficacy of online documentation is demonstrated by the growing 
number of systems that offer it in some form. Good online information 
offers convenience and timeliness to the end user and flexibility and 
ease of updating to the documenter. This study examines and evaluates 
current LaRC online documentation capabilities, with specific 
attention to: 


o NOS online Help (as part of the CDC system) 

o Information Tree (IT) — a Langley generated 
program for allowing users to keep up with 
programs and capabilities developed by the LaRC 
community 

o KEYWORD — a listing and description of utility 
programs available 

o CONTEXT — a CDC "tool" for constructing online 
documentation 

In addition, the study summarizes state of the art by reviewing online 
capabilities such as the Symbolics "Document Examiner," UNIX's "man" 
and "apropos" features, and the ENACS text editor's INFO system. 

Summary of Recommendat i ons 

The industry trend is toward more online documentation. Furthermore, 
many installations have begun converting offline to online. However, 
simply putting existing texts online is not a satisfactory approach 
because accessibility and ease of use are seriously impaired. Two 
general approaches help to overcome these limitations: (1) design 

and rewrite the manuals so that substance, structure, style, and 
syntax meet the requisites of online presentation; (2) treat existing 
documentation as a database and develop sophisticated interfaces 
allowing the user access to appropriate Information in a timely 
manner. This report adopts a hybrid of the two methods and suggests a 
scenario for online documentation at LaRC, phased in accordance with 
predicted hardware and software advances. 
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COMPUTING THE EXPONENTIAL FUNCTION OF A 
BI-DIAGONAL MATRIX WITH APPLICATIONS 
TO RELIABILITY ASSESSMENT OF 
FAULT-TOLERANT COMPUTER SYSTEMS 


Robert L. Causey 
Professor of Computer Science 
Christopher Newport College 
Newport News, Virginia 23606 

ABSTRACT 


The Semi-Markov Unreliability Range Evaluator (SURE) software package contains a 
module which estimates the probability E(T) of arriving at the death state of a 
pure Markov process by a given (mission) time T. The pure Markov chain is a sub- 
path of the semi-Markov model which describes the behavior of fault-tolerant 
computer architectures. 

The pure Markov model is characterized by a system of n first-order ordinary 
differential equations where n is the number of states leading from the intial 
state S^ to the death state S R . If p^(t) denotes the probability of being in 
state at time t, the differential equations are of the form p, f (t) = a^p^(t), 

P n ' (t) " iP n _i(t) and p^Ct) = ^i_ 1 P i _ 1 (t) + for 1 ■ 2, .... n-1 where 

XjL denotes the transition rate between state S. and state . The constant ot^ 

is given by -A^ -y^ where y^ denotes a stransition rate from to an intermediate 
state in the semi-Markov model. There are exactly n-1 as and As. We note also 
that a^ <.0 ; A^ > 0 for all i. 

The system of differential equations may be expressed in the matricial form P T = AP 
where P ® P(t) denotes the colum matrix P(t) = Cp-^t), ••• >P n (t)^ * (* denotes 
trabspose) and A = (a . ) denotes the n by n bi-diagonal matrix with a^ = 
a^^ i = A^ (i = 1, ,n-l) and a^ . = 0 elsewhere. The eigenvalues of A are 0, 

a^, ,a n , . The probability function E(T) is calculated by solving the matricial 

system P 1 = AP as an initial-value problem with P(0) = 1_1 , 0, ... ,CQ = e * and 
then setting E(T) = P n (T). The general solution is given by P(t) = exp(At)e. where 
exp(M) denotes the matrix exponential of M. The value of E(T) is therefore the 
(n, 1) component of exp (AT) . 

The primary objective of the summer research has been to develop and implement fast 
algorithms for computing E(T) . The following results may be reported. Let exp(AT) 
=[^Cj, ... where c^ denotes the ith column of the exponential. We have proved 

that c n = e n = [O, 0, .^. , l]* and that 

(1) (A - a i I)c i = for i = 1, 2, ... ,n-l. 

These relations hold regardless of the multiplicity of the eigenvalues of A or 
whether A is defective. The coefficient matrix in (1) is always singular of rank 
n-1. Consequently, the solution space is one-dimensional. If the eigenvalues of 
A are distinct, the equation (1) may be used to solve for the columns of exp(M) 
t^here, M = ATT} recursively in the order c ,, c ... ,Cj by noting that exp(M) 

is lower triangular with the scalar exponentials exp(a^T) along* the main diagonal 
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in the first n-1 columns. Software is being developed to implement this algorithm. 

If A has multiple eigenvalues, the algorithm requires modification. This problem 
is the subject of further investigations. 

Another module for the fast calculation of E(T) is being developed. This has 
complete generality and is based on the rational approximation to exp (AT) found in 
the ninth diagonal entry of the Pade table for exp(x). The code takes full 
advantage of the sparse form of A. While it is essentially equivalent mathematically 
to the algorithm used in the NASA subroutine CONEXP, it will be far more efficient 
in its use of computing resources. 
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CONSIDERATIONS IN A COST/BENEFIT ANALYSIS 
OF THE COMMERCIALIZATION OF SPACE PROGRAM 


Bernadette P. Chachere 
Department of Economics 
Hampton University 
Hampton, VA 23668 


Abstract 


In his state of the union message of January 1984, President Reagan announced 
the Implementation of a number of executive Initiatives designed to promote 
private sector Investment In space. The executive Initiatives have been termed 
the 'commercialization of space program. 1 ' A primary Instigation of this 
commercialization of the space program can be traced to an Industry Commercial Space 
Group Report of December 1983 which concluded "commercial activities In space by 
private enterprise need to begin now If our nation Is to retain Its leadership 
In science and technology. . .and Its advantages In international trade. With 
government as a partner, private enterprise can turn space Into an arena of 
Immense benefits for our nation." More specifically the macroeconomic goals of 
the space commercialization program as articulated by NASA Commercial Space Policy 
(unpublished data, NASA Commercialization Office, 1985) are to promote industry 
participation in accelerating the application of space technology to the civilian 
sector to "strengthen the nations economic vigor; retain our international 
technological preeminence; improve our balance of trade and maintain our high 
standard of living." 

This study assumes that the significant and positive correlations found In 
previous studlesl between technological advance, economic growth, trade 
advantages and Improvements In general economic welfare will hold to a greater 
or lesser degree In the application of space technology to civilian goods and 
services. The stud(y also assumes that technological progress can be Influenced 
directly by the resources committed to research and development activities, 
l.e. technological progress/change Is not a random event. What the study 
questions Is whether or not there will be a greater amount of resources 
committed to RAD projects In the space arena under the 'commercialization' 
program Inltlatves than there would be In the absence of the initiatives. 

Commercialization means to make profitable In a business way.? T hus the 
assertion of the commercialization of space program Is that technological 
progress will be accelerated If the allocation of resources to RAD activities 
takes place In the private, for-profit sector, rather than the public, 
non-profit sector. Thus the national Interest will be better served via the 
exploitation of space rather than the exploration of space 3 . 
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Evaluating the potential public or social benefits of the commericallzation of 
space program requires answers to the following questions: 

1. What determines the level of a firm's research and development 
expenditures? While the term 'industry' is used repeatedly In 
discussions of the commercialization program. Industry as such does not 
constitute a legal nor decision making entity. Allocation of resources 
takes place at the level of the individual firm. 

2. What determines the allocation of a firm's expenditures among 
alternative RAD projects? Alternatively stated, what are the 
characteristics of the RAD projects that are undertaken? Social 
benefits are derived not only from the quantity of R&D expenditures 
but also from their quality and relevance. 

3. What determines the rate of diffusion of technological change? Full 
social benefits of innovations are not realized unless the new process 
or product is adopted by other firms. It is the imitation by other 
firms which moves technological progress to the industry level and 
enables its impact on macroeconomic variables. 

4. To what extent can the specific initiatives of the space 
commericalization program influence positively the quantity and 
character of RAD activities of firms and industries? 

The analysis of the study uses the theoretical and econometric models developed 
in the economics literature to isolate the key determining variables relevant to 
the questions posed in question (1), (2), and (3) above. 4 Specific 
commercialization program initiatives are than hypothesized as to their 
positive, negative or neutral effect on the total RAD expenditures of individual 
firms, the allocation of expenditures among alternative RAD projects and the 
intra firm rate of diffusion of RAD outcomes. The analysis is speculative given 
the lack of data and the vague articulation of the implementation of the 
initiatives. Both undoubtedly reflect the newness of the program. 


57 



NOTES 


1. This relationship has been documented In numerous studies employing various 

methodologies. The more definitive Including the following: E. Denison, 

The Sources of Economic Growth in the United States and the Alternatives 
Before Us (New York: Committee for Economic Development, 1962); R. Solow, 
’'Technical Change and the Aggregate Production Function," Review of Economic 
and Statistics, August 1957, Technology and the American Economy, Report of 
the National Commission on Technology 7 Automation and Economic Progress 
(Washington, 1966). 

2. Webster's New Collegiate Dictionary 

3. Exploration Is defined as a systematic search for discovery. Exploitation 
is defined as utilization or to obtain value for one's own advantage or 
profit. Exploitation Is the word used by the National Space Policy 
published on July 4, 1982 

4. The models used draw heavily on the seminal work of E. Mansfield, Industrial 

Reseach and Technological Innovation (New York. Norton & Co., 1968) and K. 
Arrow, "Economic Welfare and the Allocation of Resources for Invention," in 
the National Bureau of Economic Research Conference Report, The Rate and 
Direction of Inventive Activity (Princeton: Princeton University Press, 

mi; 
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A STUDY OF POTENTIAL GASES FOR 
BLACKBODY-PUMPED TRANSFER LASER 


Kwan-Yu Chen 
Department of Astronomy 
University of Florida 
Gainesville, Florida 32611 

Abstract 

The conversion of solar radiation into high- power laser 
radiation in space could have wide application in future 
space power and propulsion requirements. One concept is to 
collect and converge solar radiation in space into a black- 
body cavity heating it to a temperature of, say, 2000 K. A 
gaseous material passing through a transparent chamber with- 
in the blackbody cavity is excited radiatively. The excited 
molecules, Mj.are then mixed with another gaseous material 
of molecules, M£. In molecular collisions, M} transfer their 
excited energy to Mo; and this creates a population inversion 
in some energy level of M 2 • This study deals with the inves- 
tigation of potential materials that might be used in such a 
transfer laser. 

The radiative temperature of the blackbody cavity would 
be in the range from 1500 K to 2500 K. The corresponding 
wavelengths of the maximum radiative intensities are 1.93 ym 
and 1.16 ym respectively. In the irfrared region of the 
spectrum, the molecules would be excited to higher vibratio- 
nal levels. The i ntermol ecul ar processes involve both vi- 
bration-vibration transfer and vibration-translation trans- 
fer of energy. 

Diatomic heteronulear molecules are considered only for 
Mj, while triatomic and other polyatomic molecules are con- 
sidered mainly for M2- A prelimary list of 12 diatomic, 10 
triatomic, 11 four-atomic, 10 five-atomic, 3 six-atomic, and 
1 seven-atomic molecules is compiled. Their fundamental 
vibrational levels are noted. Large rate constants for vi- 
bration-vibration transfers between Mj and M 2 a re desired 
for higher laser efficiencies. In these cases, their vi- 
brational levels should be matched closely in order for 
efficient resonant transfer to be acheived. 

To illustrate the principles of the mechanics of the 
transfer laser operation, an example of laser-excited infra- 
red laser with Mj=C0 and M2=C02 is discussed. 


59 



THE ANALYSIS OF HYDROCARBONS 
IN THE TROPOSPHERE 


Randolph A. Coleman 
Prof BMor 

Department of Chemistry 
College of William and Nary 
Williamsburg, Virginia 231SS 

An analytical system has been developed for the 
determination of trace volatile organic compounds in the 
atmosphere. Current testing is focused on the analysis of 
hydrocarbons in the troposphere as part of NASA's efforts to 
gain an understanding of tropospheric chemistry. As 
presently configured, the system utilizes a per easel set ive 
membrane (Dupont Naf ion) for water removal, a low-pressure 
liquid nitrogen cryogenic trap for concentrating the sample, 
and a gas chromatograph with a flame ionization detector for 
separation and quantification of components in the sample. 
Bas sample sizes of 300 oL are sufficient for the 
determination of components at concentration levels as low 
as sub-part— per — bi 1 1 ion volume mixing ratios. 

To date the system has been used successfully in the field 
for the analysis of hydrocarbons in samples taken from 
40,000 ft down to sea level. Altitude/concentration profiles 
were determined for normal alkanes and lower molecular 
weight aromatics. In addition to these atmospheric samples, 
the system has also been used to analyze trace hydrocarbon 
co mp one n ts in the nearly pure methane gas emitted from the 
anaerobic sediments found in freshwater swamps, salt 
marshes, rivers, and coastal waters. The reduced-pressure 
cryogenic trapping loop in the system proved to be very 
effective for the latter analyses, since methane is not 
concentrated in the loop, but higher molecular weight 
components are. Continued testing this summer will also 
include the analysis of soil gas and leaf litter 
decomposition products from forest floors. 



TRANSONIC FLUTTER ANALYSIS 
BY A LEAST-SQUARES FINITE ELEMENT SCHEME 


Christopher L. Cox 
Professor 

Department of Mathematical Sciences 
Clemson University 
Clemson, SC 29631 

Abstract 

The problem under consideration is unsteady transonic flow about a harmonically 
oscillating airfoil, where the unsteadiness is taken to be a linear perturbation 
about a nonlinear steady state. This assumption plus harmonic time dependence 
results in a linear equation for a complex potential <|> 

(|_1^m2 - (y + 1) <i>o x J4’x)x + fyy - 
M 2 t 

Zi^ 1 4>x + 1 ^ ~ 1“ (y - 1 ) 't’Oxx J + = 0 (1) 

e e 

having discontinuous coefficients which depend on the solution of the steady state 
equation 

[JL^M? - (y + 1) <j>o x J 4>0 X X + ’•’Oyy = 0 (2) 

M 2 e 

In these equations M is the free stream Mach number, c is a scaling parameter 
depending on M and airfoil thickness, y is the ratio of specific heats and u> is the 
reduced frequency. 

The test airfoil is a 6 % thick parabolic arc and we solve for M = .872 and 
pitching amplitude of .1 degree. Symmetry considerations allow us to solve in the 
half plane shown in figure (1). In addition to equation (1), we must impose 
appropriate boundary conditions along the airfoil surface, wake, and outer 
boundary. The steady state potential 4>o is computed using the finite difference 
code TSFOIL. 

Examination of the coefficient term 

CF - J-H2 - (T ♦ 1) *„ x (3) 

Vr e 

provides an understanding of the main characteristics of this flow. Figure (2) is a 
plot of CF along the airfoil surface. When the curve first crosses the x-axis 
equations (1) and (2) change type from elliptic to hyperbolic and the flow changes 
from subsonic to supersonic. These characteristics are reversed abruptly near x = 
.75 signifying the root location of a shock, a line in the flow plane along which a 
discontinuous jump occurs in 4> 0 and CF. 
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A least squares approach is convenient for this problem because there is no type 
of dependence like that in finite difference methods, and boundary conditions are 
imposed with ease. In addition, the finite element grid is quite flexible and points 
can be clustered locally in regions of large gradients while leaving the grid coarse 
in the far field. A typical grid is shown in figure (3). Points are clustered in 
the two regions where singular behavior occurs, around the leading edge at (0, 0) and 
near the root of the shock at (.75,0.). 

Since the mesh used in TSFOIL consists solely of horizontal and vertical lines, 
the fitting of <j> 0 to the finite element grid points demand particular attention, 
especially near tne shock. The difference scheme spreads the shock somewhat over 3 
or more points whereas a discontinuous jump between two successive points is more 
reasonable. Thus it is necessary to fit the shock by extrapolating from either side 
so that there is a clean jump in <j> 0 . Figure (4) displays the fitted CF versus 
the finite difference data, at y = u. A subroutine was written this summer to fit 
(|> 0 and <|> 0 at each row of points which passes through the shock in a simple 
manner. tR£ interpolation scheme for fitting <j> 0 and <j> 0 to the finite 
element grid was revised to ensure the jump at tne shock x is preserved. 

Because the derivatives of <f> are of as much or more physical significance that $ 
itself, we solve a first order formulation of (1). The formulation used previously 
to this summer is 

u x + Vy - 2i(o <j> x + (o£_ - i u) ( y - 1 ) <t> 0xx ) <p = 0 ( 4 ) 

£ £ 

u - CF 4> x = 0 

v - <py = 0 

This choice for the variables was deemed appropriate as it assured necessary 
conti nuity- in the product of two discontinuous variables through the shock. Of 
primary interest is the jump across the airfoil in the pressure coefficient, which in 
linearized form is 


Cp = -2 (<f>x + 1H>) 


(5) 


The pressure jump computed using equation (4) is shown in figure (5). A point 
of controversy among those studying this problem is whether or not the formulation 
accounts for shock notion, which is assumed to be harmonic. The apparent answer is 
no. Thus an additional condition must be imposed, namely 

< 2iu> + CF X > [<|>J - [CF J <<j) X > = 0 (6) 

£ 

where < . > denotes average value and [ . J denotes the jump across the shock. This 
necessitates building the shock into the finite element grid as shown in figure (6). 
It is appropriate then to let the variable u represent <p x since equation (6) 
contains the necessary continuity conditions. 

Computer experiments are now being performed with this formulation to determine 
the optimal size of computational region and number of grid points. No analytical 
solution exists so comparisons will be made with experimental results and results 
using other methods. 
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SURVEY OF THE QUALITY OF THE 
RESEARCH ENVIRONMENT AT NASA LANGLEY 


Anthony Dalessio 
Department of Psychology 
Old Dominion University 
Norfolk, Virginia 23508 


Background 

The purpose of the project was to assess the quality of the research 
environment at Langley Research Center (LaRC). Specifically, the project was 
an effort to identify the aspects of the research environment which the LaRC 
scientists and engineers find to be of high quality, and the areas which are 
in need of improvement. Feedback on these topics will be provided to the 
Langley management. To more clearly define the aspects of the research 
environment on which to focus the study, interviews were first conducted with 
members of the LaRC senior staff. The areas of interest defined in these 
conversations were: (1) the amount of challenge, variety and freedom in the 
work assignments of the scientists and engineers; (2) the quality of the tools 
and support; (3) paperwork; (4) the research guidance process; (5) the NASA 
publication process; (6) communication within and between Langley 
organizations; (7) the ability of the scientists and engineers to initiate 
projects; (8) credit, recognition, and visibility of work; (9) career 
planning; (10) job satisfaction. 


Study Design 

Two methods were used to obtain data in these areas. The first method was a 
questionnaire distributed to 312 scientists and engineers in 5 LaRC 
Directorates (Electronics, Structures, Aeronautics, SE&O, and Space). 
Questionnaire items dealing with challenge, variety, freedom, job satisfac- 
tion, supervision, promotions, and awards were drawn from previously developed 
standardized questionnaires, and from questionnaires previously distributed by 
the Office of Personnel Management and NASA headquarters to LaRC scientists 
and engineers. Questionnaire items in other areas were generated by the 
author for the purposes of this study. Currently, completed questionnaires 
have been received from 184 of the employees. Data are presently being coded. 
Analyses will include presenting overall Center means and Directorate means 
for responses to the items. Comparison of present responses with those of 
past responses to the Office of Personnel Management and NASA headquarters 
items will also be made. 

The second method used to collect data was in-depth structured interviews. 

The opinions of fewer researchers could be examined with this method compared 
to the questionnaire method; however, more detailed information could be 
gathered from each respondent. Interview questions were developed to assess 
the areas previously noted. A total of 23 researchers representing the five 
Directorates have presently been interviewed. All participants were asked the 
same set of questions and to elaborate on responses where approporiate. 

Each interview lasted approximately one hour. Interview participants 
were randomly selected from the pool of people who responded to the 
questionnaire. 
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Preliminary Results and Conclusions 

Based on the interview data gathered thus far, several preliminary conclusions 
can be drawn about the research environment at LaRC. 

(1) The scientists and engineers generally view their jobs as providing them 
with the opportunity to test new innovative ideas and to develop new methodo- 
logies. A number of respondents cited these opportunities as the major 
advantage of the job. 

(2) The paperwork associated with submitting manuscripts for outside review 
at journals and conferences was viewed as redundant and unnecessary in many 
cases. Also the researchers felt too much time was spent in justifying the 
research that they were doing through paperwork and presentations. 

(3) There was some mixed reaction to the tools and support available for 
conducting research. The LaRC Technical Library was described as excellent by 
a number of the respondents. Computer support was seen as adequate by some 
and inadequate by others. Generally, employees who had computers in their own 
branch were satisfied with their systems. The central computing facility was 
viewed as adequate, but turnaround time! was described as problematic. 

(4) Most of the scientists and engineers felt they had adequate opportunity 
to initiate research projects provided funding was available and the project 
was consistent with the goals of their branch. Ability to initiate projects 
relevant to work being conducted in other branches was thwarted. The 
scientists and engineers felt they had complete freedom to choose the methods 
they used to address their research problems. 

(5) The participants felt that they were undercompensated relative to 
industry. Most were willing to accept the compensation differential for the 
opportunity to participate in the unique work at LaRC. However, many felt 
that the compensation gap between industry and NASA was growing larger, and if 
the gap became too large or if the research environment changed, the younger 
scientists and engineers may begin to look for opportunities in other 
organizations. 

(6) Many of the scientists and engineers felt that the reviews provided by 
NASA technical review committees resulted in high quality NASA Technical 
Papers. There was not complete agreement on this issue however. Several of 
the respondents felt that at times papers were not given a thorough review 
through the NASA publication process. This situation was seen as prevalent 
when reviewers were chosen by the researcher, or when reviewers did not have 
the expertise to review the paper. 

(7) The scientists and engineers generally felt that there was a good amount 
of technical communication within the branch. They also felt that there was 
not enough technical Information shared between branches working on similar 
problems. Finally, the researchers felt that greater support for travel would 
greatly Increase their contacts with professionals outside of LaRC. 
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MAKING SENSE OF FLIGHT TRANSITION MEASUREMENTS 


Henry P. Day 
Assistant Professor 
Engineering Science 
Virginia Polytechnic Institute 
and State University 
Blacksburg, VA 24061 


In a recently completed flight experiment, a Cessna Citation was used to 
measure boundary-layer transition in flight. The present report discusses 
that experiment and makes some recommendations for an upcoming NASA experiment 
using a Lear Jet. The purpose of the Citation experiment was to find how much 
Tollmien-Schlichting waves and cross-flow vortices interact in causing transi- 
tion. The results were disappointing and do not allow any conclusions. The 
present author did several analyses to improve the results or explain the 
scatter. The effect of pressure measuring belts on pressure measurements was 
studied briefly, though the proper computational tools are not available. The 
airfoil and boundary-layer analysis codes available at Langley were run to 
compare with Boeing codes, in particular to determine if laminar separation 
could be blamed for transition in some cases. Finally, and most usefully, the 
difference between moving and stationary cross-flow vortices was checked to 
see if their motion must be accounted for. Better cross-flow vortex codes 
than are presently available will be needed to adequately interpret flight 
data. 
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COMPUTATION OF TRANSONIC FLOW OVER 
BODIES IN A WIND TUNNEL 


Michael 

Associate 

Mechanical 

Valparaiso 

Valparaiso, 


L. Dona 
Professor 
Engineering 
University 
Indiana ^6383 


One of the goals of computational fluid dynamics is to supplement the wind 
tunnel as an aerodynamic design tool. Therefore it is desirable to have 
computational models for flow over bodies in a tunnel. Such models are 
useful for determining corrections to wind tunnel data. 

During my previous experience at NASA Langley I developed a grid generation 
program for two dimensional bodies of general shape in a tunnel with shaped 
walls. The program uses a series of Schwarz - Christoffel transformations 
and shearing transformations to map the physical region between the body and 
tunnel walls into a rectangular computational domain. This procedure leads 
to an "O-type" mesh which is nearly orthogonal everywhere. The mesh is said 
to be O-type because one set of coordinate lines wraps around the body. The 
other set originates on the body surface and terminates on the tunnel walls. 
The resulting grid has desirable properties for the numerical solution of 
the inviscid flow equations. 

The goal of my work this summer was to incorporate the tunnel grid coor- 
dinates into an existing Euler equation flow solver and hence produce a 
computational model for transonic inviscid flow over two dimensional bodies 
in a tunnel. 

The Euler code chosen was FL052, which was developed by Antony 
Jameson ( Ref erence 1) and vectorized by J. Thomas and M.D. Salas of NASA 
Langley. The code solves the unsteady Euler equations using a finite volume 
scheme for spatial differences and a multistage Runge - Kutta scheme for the 
time step. Multigrid is used to advance the solution rapidly to steady 
state. 

Early on in the work it was discovered that the grid generation program was 
not working completely satisfactorily. When it was applied to airfoils with 
blunt leading edges the scheme tended to produce small wiggles in the coor- 
dinate lines near the nose region. It was found that this was due to 
difficulties in the spline routine used to fit a smooth curve to the airfoil 
shape from the discrete input coordinates. A convenient way to fix this 
problem is to change the parameterization of the airfoil coordinates from y 
versus x to x versus s and y versus s, where x is distance along the chord, 
y is distance normal to chord, and s is arclength along the surface. The 
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arclength parameterization avoids the singularity in derivative that occurs 
in y versus x at the leading edge. The spline routine had no difficulty 
fitting smooth curves to the two well behaved functions x(s) and y(s). The 
grid generation routine now works quite well and can treat airfoils of very 
general shape . 

To incorporate a wind tunnel geometry into an Euler code written with free 
air boundary conditions in the far field, several considerations have to be 
taken i i nto account. First a solid wall boundary condition replaces the free 
air boundary condition at the outer boundary. Inflow and outflow boundary 
conditions are needed at upstream and downstream infinity in the channel. 
Finally the calculation of the cell volumes has to be modified for cells at 
upstream and downstream infinity. 

At this writing it is too early for any results but it is hoped that the 
resulting code will prove to be a useful tool for modeling transonic flow 
over bodies in a wind tunnel. 


Reference 1: Jameson, A., "Solution of tne Euler Equations by a Multigrid 
” Method , " Applied Mathematics and Computation , Vol .13,1 983, 

pp. 327-356. 



EQUIVALENT CONTINUUM ANALYSIS OF PROPOSED LATTICE STRUCTURES 

FOR THE SPACE STATION 


John 0. Dow 
Associate Professor 
Department of Civil, Environmental 
and Architectural Engineering 
University of Colorado 
Boulder, Colorado 80309 


This research develops a procedure for creating an 

equivalent continuum beam model for lattices of the type being 
proposed for the space station. The resulting equivalent beam 
models allow the following areas to be studied: 

1. Stiffness comparisons of various configurations. 

2. The relative coupling between the deformation modes. 

3. The effect of joint stiffness on the stiffness and coupling 
parameters . 

4. The effect of bi-lmear joint stiffness on stiffness and 
coupling parameters. 

5. The effect of non-uniform degradation or structural damage 
on stiffness and coupling parameters. 


The equivalent continuum analysis is based on creating an 
equality between the strain energies of the lattice structure and 
the desired continuum model. The theoretical basis of this 
research is similar to the approach developed by Noor, Anderson, 
and Greene which utilizes an a-prion defined series expansion 
relating displacements and strain gradient terms (1) . This work 
refines the approach developed by Dow, Feng, and Bodley which 
allows methods from linear algebra to be used to identify a set 
of linearly independent strain states from a candidate set in the 
displacement-strain gradient expansion (2) . The development 
presented here is specialized to certain proposed space station 
configurations. This specialization to specific structures 
allows the method of Ref. 2 to be utilized in a more physically 
based and intuitive manner. In other words, the set of linearly 
independent strain gradient terms can be identified by 
inspection. The resulting compactness of the procedure makes it 
particularly suited for the preliminary design stage and 
implementation on personal computers. 
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The relationship between the nodal displacements and nodal 
rotations as a function of strain gradient quantities are 
developed for the two dimensional case in the first step. The 
two dimensional case contains enough detail to allow the results 
of the three dimensional development to be understood. The three 
dimensional results are included without proof. The coefficients 
of the displacements, rotations, and strains are given Tables 1- 
3. An interesting by-product of this derivation is the 
generation of the compatibility equations from elasticity. This 
can be seen in Table 3. 

The second step identifies the linearly independent strain 
gradient quantities that characterize the space station lattice 
structure. Both, the case with hinged and the case with fixed 
joints, are considered. The terms identified by inspection are 
numerically verified for a candidate space station lattice. 

The strain gradient parameters which characterize a beam 
with bending around two axes, rotation around the longitudinal 
axis, two orthogonal shear deformations, and longitudinal 
extension are then defined. These variables are called kinematic 
variables and they allow the first order stiffness and coupling 
parameters to be extracted from the stiffness matrix of the 
lattice structure. The remaining linearly independent variables 
can be included in the analysis to add the second order effects 
to the stiffness and coupling parameters. 

With the identification of the linearly independent strain 
gradient quantities associated with the proposed space station 
lattice and the parameters that characterize a beam, the 
equivalent continuum beam parameters can be computed. This is 
done by transforming the strain energy in terms of the finite 
element stiffness matrix and nodal quantities to an expression in 
terms of the strain gradient quantities using the relationships 
developed in step 1. If only the kinematic variables are 
included, the first order stiffness and coupling parameters 
result. The secondary effects can be added in the beam 
parameters by using a static condensation procedure to implicitly 
include the secondary strain gradient parameters in the analysis. 

The method is applied to a candidate structure and the 
effects of the second order terms are computed and discussed. 
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TABLE 1 THE SIXTY COEFFICIENTS FOR DISPLACEMENT FUNCTIONS 



a 1 for u(x,y,z) 


q ♦ y x2 /2 
e x,x /Z 


^ Y xy,y * e y,x^ 2 
^ Y xy,z * Y yz,x * Y xz,y^ Z 
^ Y xz,z * 'z.x> /2 


e x.xy /2 

'x.xz' 2 

^ y xy,yy ~ c y,xy ^ 6 

€ x,yy fZ 

* Y xy,yz ~ *y,xz^ 2 
^ Y xz,zz " *z.xz>'« 




b,for v(x,y,z) 


V 2 + r 


(T xy.x * c x,y )/2 


^ Y xy,Z + Y yz,x ' Y xz,y^ 2 

£ y.y /z 


^ T yz,z ‘ e z,y) /Z 
^ Y xy,xx ’ e x,xy^ 6 


^ Y xy,xz * e x,yz^ Z 

c y.yy ' 6 

e y.*y /2 

e y,yz /Z 

^ Y yz,zz " e z,yz^® 

(y - e )/2 

yz.xz z,xy 


c i for w{x,y,z) 


Y xz ' 2 - « 


Y yz /2 + P 


^ T xz,x * e x.z >' 2 
^~ Y xy,z + Y yz,x + 


(T yz.y - E y.z ,/2 


e z ,z ^ 2 

^ Y xz,xx ’ £ x.xz >' 6 
Wz.xy * £ x,yz )/Z 
£ z,xx /Z 

^yz.yy * V.y *^ 6 

* Y yz,xy ’ V.xz^ 2 

W 2 


Y xz,y ^ 2 


e z,zz /6 

£ z,xz /2 

'z.yz ^ 2 


u(x,y,z) “ *j ♦ «2* ♦ *3X * * 4 * + a 5 X 2 ♦ 
w(x,y,z) • c, ♦ CjX ♦ Cjy ♦ c,,z ♦ c 5 x 2 + 


♦ * 2 Qxyz ; v(x,y,z) ■ bj ♦ b 2 x + bjy + b„z ♦ b s x 2 ♦ ... + b 20 xyz j 
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TABLE 2. THE COEFFICIENTS FOR ROTATION FUNCTIONS 


COL. 

# 

TERM 

(1) 

p(x.y.z) 

(2) 

q(x.y.z) 

(3) 

r(x,y,z) 

(4) 

1 

1 

P 

q 

r 

2 

X 

(y “Y )/2 

VT xz,y T xy,z" 

<-W 2 ‘x.z>' 2 

(l xy.*' 2c x.y )/2 

3 

y 

U yz,y- Ze y,z )/Z 

^ Y xy ,z“ Y yz,x^ 2 

( ' r xy.y t2i y.x )/2 

4 

Z 

(" Y yz,z +2e z,y) /2 

<W 2t *.« )/2 

^~ Y yz,x~ Y xz,y^ 2 

5 

X 2 

( v -v ) /4 

v 'xz ,xy 'xy .xz " 

^ ~ Y xz ,xx +2e x,xz )/4 

^ Y xy ,xx" 2e x,xy ^ 4 

6 

xy 

^ Y yz,xy~ 2e y,xz^ 2 

^” Y xz,xy +2e x t yz^ 2 

( e - e )/2 

y.xx x,yy" 

7 

xz 

^ _Y yz,xz +2e z,xy^ 2 

^ e x,zz" e z,xx^ 2 

( V.xz' 2e x,yz )/2 

8 

y 2 

^ T yz,yy' 2e y,yz^ /4 

^ Y xy.yz” Y yz,xy^ 4 

^’ Y xy»yy t2 ‘V.yx^ 4 

9 

yz 

^ e z,yy" e y,zz^ /2 

^ Y xz ,yz” 2e z,xy^ 2 

( - Y +2 e ) /2 

T xy,yz y,xz'' 

10 

Z 2 

f-v +2e ) /4 

v 'yz.zz z,yz" 

^ Y xz,zz“ 2e z,xz^ 4 

( y - y )/4 

VT yz,xz T xz,yz" 




TABLE 3. THE COEFFICIENTS FOR STRAIN FUNCTIONS 


TERM 

e x (x,y.z) 

e y (x.y.z) 

e z (x,y,z) 

Y xy (x,y,z) 

Y yz (x,y,z) 

Y xz (x.y,z) 

1 

E x 

E y 

E z 

Y 

T xy 

T yz 

Y xz 

X 

e x,x 

€ y.x 

e z,x 

Y 

xy ,x 

Y yz,x 

Y xz,x 

y 

e x,y 

e y.y 

e z,y 

T xy.y 

Y yz,y 

Y xz ,y 

z 

e x,z 

e y.z 

E Z.Z 

Y xy,z 

Y yz,z 

Y xz,z 

X 2 

l/2e 

x,xx 

1/2e y.xx 

1 / 2 e z,xx 

l/2Y xy,xx 

»— * 
ro 
** 

X 

1/2 y 
' xz,xx 

xy 

€ x,xy 

E y.xy 

E z,xy 

Y * 

xy , xy 

Y yz,xy 

Y xz,xy 

XI 

€ x # xz 

E y.xz 

E z ,xz 

Y xy ,xz 

Y yz ,xz 

Y * 

T xz ,xz 

y 2 

i/2 c 

x.yy 

l/2e 

y.yy 

1/2e 

z,yy 

l/2y 

xy.yy 

1/2 y 

xy ,yy 

1/2y * 

xz,yy 

yz 

€ x,yz 

E y.yz 

E z.yz 

Y xy,yz 

Y yz, yz* 

Y xz ,yz 

z 2 

l ' 2 ‘x.zz 

1/2E y,zz._ 

1 ' 2 e z,zz 


1/2Y yz.zz 

^xz.zz 


* The terms starred are equal to the compatibility relations which are given below. 


Y 3 

xy ,xy 

e x,yy 

♦ 

E y .xx 

Y xy.zz 

* Y 

xz ,yz 


Y yz,xz 

’ 2e z,xy 

Y yz,yz ' 

e y.zz 

♦ 

E z.yy 

Y yz,xx 

■ Y 

f xy,xz 

+ 

Y xz ,xy 

2e x,yz 

Y 3 

xz.xz 

€ x,zz 

+ 

f z ,xx 

Y xz ,yy 

* Y 

xy *yz 

+ 

Y yz,xy 

' 2e y.xz 
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MECHANICAL AND EPR* STUDIES OF THE RADIATION-DURABILITY 
OF FIBERS AND FIBER/RESIN COMPOSITES 


Milton W. Ferguson 
Assistant Professor of Physics 
Department of Chemistry /Physics 
Norfolk State University 
Norfolk, Virginia 23504 


The use of fiber-reinforced composite materials for space structural 
components has greatly increased in recent years. The principal advantages 
of these materials are high strength, high modulus of elasticity, and low 
density. Another advantage is the ability to tailor the materials to meet 
specific load and stiffness requirements. However, in aerospace applications, 
these composites are exposed to solar ultraviolet radiation and charged 
particles (protons and electrons). Therefore, if these composites are to 
be used as efficient high performance materials, their long term radiation- 
durability must be established. 

The performance of fiber-reinforced composite materials depends not 
only on the strength of the component materials but also on the interaction at 
the fiber resin interface. The interfacial bonds, which influence the ultimate 
stress, modes of failure, interlaminar shear strength, and modulus of elas- 
ticity of the composites, are sensitive to radiation. Exposure of the 
structure to electron bombardment induces free radicals which can lead to 
profound changes in these mechanical properties of the material. 

The goals of the study were: (l) to determine and analyze the effects 

of radiation on the mechanical properties of graphite and Kevlar fibers; 

(2) to determine and analyze the effects of radiation on the mechanical 
properties of graphite/resin composites fabricated from the same types of 
graphite fibers; (3) to see if a correlation exists between the radiation 
durability of the component materials and the fiber/resin composites. 

Graphite (Celanese Celion 6000 sized with 0 . 1 % polyimide, Celanese Celion 
6000 nonsizgd) and Kevlar 49 fibers were exposed to 1 MeV electrons up 
to 9.5 X 10° rad. Since the results of investigations on the radiation- 
durability of Ultem, a polyetherimide , have been previously reported by my 
NASA associates, Drs. Edward and Sheila Long, Celion 6000 nonsized graphite/ 
Ultem and Celion 6000 sized graphite/Ultem composites were fabricated. 


♦Electron Paramagnetic Resonance 
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These composites, with fibers oriented uniaxially in the longitudinal 
and transvers directions, were exposed to the same dosage of radiation 
as the isolated fibers. Characterization of the fibers and composites 
was performed by tensile tests and interlaminar shear strength tests, 
respectively, before and after irradiation. Electron paramagnetic 
resonance (EPR) data were used to detect, identify, and determine 
the densities of radicals created by the irradiation. 

The mechanical data indicate that there were radiation induced 
changes in the tensile properties of the graphite and Kevlar 49 
fibers. These results are supported by the presence of induced free 
radicals as shown by the EPR measurements. Interlaminar shear strength 
tests and EPR analysis of the composites are continuing. 
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SAMPLING AND ANALYSIS OF 
STRATOSPHERIC AEROSOLS 


Thomas A. Gosink 

Research Associate of Geochemistry 
Geophysical Institute 
University of Alaska 
Fairbanks, Alaska 99701 

ABSTRACT 

A strategy for sampling and analysis of stratospheric aerosols will be 
discussed. The objective is to sort and strip out the various natural 
(oceanic, crustal, meteoric) and anthropogenic pollution components so that 
volcanic injections to the stratosphere can be determined separately. The 
ultimate goals are to follow the transport, transformation, and removal of 
these materials and to evaluate their climatic effect. 

Rapid methods of analysis, such as x-ray emission spectroscopy coupled 
with scanning electron microscopy, have been commonly employed in the past. 
These methods are quite acceptable for rapid first approximations, but for 
the above objectives and goals, much more detailed chemical analyses are 
required. 

Recommendations include obtaining larger samples, separate samples for 
chemical analyses and analyses using more sensitive methods of analyses 
whenever possible. 



THE INSTRUMENTATION OF A MILLIMETER 
WAVE COMPACT RANGE 


Frank B. Gross, III 
Assistant Professor 
Electrical Engineering Department 
Florida State University 
Tallahassee, Florida 32306 


ABSTRACT 


The basic principle behind compact ranges is discussed along with the 
attendant complexities and problems. A millimeter wave system is proposed 
that will help in reducing unwanted clutter and thus will better simulate 
outdoor conditions. The entire system is then controlled by a PDP11/23 
computer in order to automate measurements and to provide on site signal 
processing . 


INTRODUCTION 


The compact range is an indoor chamber dedicated to measuring the radiation 
pattern of antennas or the backscatter from radar targets. Millimeter 
electromagnetic waves range from 30-300 GHz. The distinct advantages of a 
compact range are that measurements can be made regardless of weather 
conditions and that measurements can be made on scale models due to the 
theorem of similitude. 

The chamber consists of the antenna feed system located at the focus of a 
parabolic reflector, the target located on a pedestal and the associated 
electronics and controller. The system is depicted in figure 1. 



Figure 1. SIDE VIEW OF COMPACT RANGE 
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The inherent problem arising from use of an indoor range is that energy 
scatters not only from the target but also from the walls, ceiling, floor 
and reflector edge. The unwanted scattering mechanisms are called clutter, 
and due to their spatial location arrive at the receiver at different times. 


MILLIMETER WAVE SYSTEM 


One method of reducing the clutter is by pulsing the radar signal 
transmitted, waiting for the pulse to travel round trip to the target and back, 
and opening the receiver (gating) at the predicted instant in time associated 
with the target location. The gate is then closed to prevent clutter, 
traveling from further sources, from entering the receiver. This technique 
is called a pulsed-CW System. The total round trip delay is: 

* ftt+JJLm 73 nsec (i) 


The source frequencies are 27-40 GHz. The pulse repetition rate is 5 MHz 
and the pulse width is 50 nsec. The pulsed-CW System is shown in figure 2. 
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THE TRANSITION FROM NACA TO NASA AT LANGLEY RESEARCH CENTER, 

A HISTORICAL PERSPECTIVE 


James Hansen 
Professor 

Department of History 
University of Maine 
Orono, Maine 04473 


Summar / 


The original parent organization of Langley Research Center was 
the National Advisory Committee for Aeronautics. Congress established 
the NACA in 1915 "to supervise and direct the scientific study of the 
problems of flight with a view to their practical solution." The NACA 
performed this assignment with distinction until 1 October 1958 — when 
a new federal agency, the National Aeronautics and Space 
Administration, absorbed and continued the aeronautical duties of the 
NACA and took on the additional responsibility of exploring outer 
space. The crucial event in the transition from NACA to NASA 
occurred on 4 October 1957s the Soviet Union shocked the world by 
launching the first man-made satellite into Earth orbit. Huddling 
near radios and televisions, Americans heard the rebroadcast "beep- 
beep-beep" of Sputnik, and they were alarmed. The distant Russian 
satellite seemed to be signalling a new communist threat to their 
national security and to be mocking the United States and its 
preeminent scientific and technological establishment for failing to 
be first into space. Less than one month later, Sputnik II repeated 
the insult, and the space race was on. 

This presentation outlines preliminary results of a historical 
investigation into the question, how did NACA Langley react to the 
Sputnik crisis'? Findings will be used to support the writing of the 
first chapter of a sequel to my book Engineer i n Charoe : A H i story of 
the Lanol ey Aeronau 1 1 cal Laboratory . 1917-1958 (in press, NASA, 
Dash i ngton , D.C.). (See Attachment for chapter outline of Engineer in 
Charoe . plus chapter outline of projected sequel.) 


Anal >• = i = 


In retrospect, it is surprising to find NACA Langley so unmoved 
by Sputnik. At the annual meeting of the Main Committee of the NACA 
on 10 October 1957 (only 6 days after the launching of the Russian 
satellite), the subject never came up. Nor was there any mention at 
all of Sputnik in Langley's in-house newspaper, the A i r Scoop . until 
the 3 January 1958 issue. The NACA reaction was not unprecedented, 
though. UJhen asked at a press conference on 9 October 1958 if Sputnik 
worried him, Eisenhower had replied, "so far as the satellite itself 
is concerned, that does not raise my apprehensions, not one iota." 
Eleven days later a member of the president's staff discounted the 
Russian achievement as an attempt to draw the U.S. into "an outer 
space basketball game." 
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Like Eisenhower, NACA leaders did not at first see Sputnik as a 
reason to upset the progress of the country's existing aircraft and 
missile programs. At a meeting of the House Select Committee on 
Aeronautics and Space Exploration in 1958 Hugh Dryden, the NACA's 
director for research, was asked if a crash program should be 
instituted to launch a man straight up into space and return him by 
parachute. (Werner Von Braun had suggested this earlier to the same 
committee as something that should be done just for the sake of 
beating the Russians.) But Dryden responded by calling the idea a 
mere "stunt" comparable to shooting a woman out of a cannon. 

The predominant attitude of the NACA during the late 1940s and 
1950s was to avoid "Buck Rogers" stuff. In 1947 a newspaper poll 
indicated that 29% of the American public expected i n terp 1 ane tary 
travel by 1960. Upon seeing the poll, Jerome Hunsaker, the NACA 
chairman, said to a colleague that all the poll really indicated was 
"that people who know about rockets like them. The same is no doubt 
true with regard to alcoholic drinks and chamber music." One 
prominent individual at Langley who was especially critical of "space 
cadets" was John Stack, chief of the Compr ess i b i 1 i ty Research 
Division. Stack resisted the idea of an expensive space program long 
after Sputnik — in the early 1960s he told some colleagues that he did 
not buy all that " to-the-Moon-by-noon " garbage. (In June 1962 Stack 
retired from NASA to become vice-president of Republic Aviation 
Corporation — a post at which he could continue to work on aeronautical 
projects.) 

Those NACA researchers who, unlike Stack, had been advocating 
increased work in astronautics were frustrated by the news of Sputnik, 
but they were perhaps even more frustrated by the overreaction of the 
American public to that news. The NACA had initiated studies "of the 
problems associated with unmanned and manned flight at altitudes from 
50 miles up, and at speeds from Mach number 10 to the velocity of 
escape from the earth's gravity" as early as 1952. Subsequently, 
researchers at Langley and the NACA's other field centers had 
concentrated increasingly upon the solution of problems involved in 
the development of missiles and hypersonic aircraft. (At that time, 
most NACA researchers considered "space flight," not as the ballistic- 
type flight to which we became accustomed in the 1960s, but as a 
natural extension of aerodynamic flight through the atmosphere into 
space and return.) But until the launching of the Sputniks, “space" 
had been a dirty word, especially in American political arenas. A 
high-ranking NACA official recalls that the NACA in the pre-Sputnik 
period stood "as much chance of injecting itself into space activities 
in a real way as an icicle had tof surviving] in a rocket combustion 
chamber . " 

Langley's first official reaction was to allay public fears by 
putting Sputnik in proper perspective. Spokesmen for the research 
center informed the local news media that, although the Russian 
satellite program posed a serious challenge to the United States, the 
U.S. was well prepared to meet the challenge. Langley, specifically, 
was already engaged in research applicable to the problems of space 
flight. Given needed additional manpower and research equipment, the 
lab would soon be doing its part to conquer space. 
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ATTACHMENT 


Outl i ne of Eno i neer i n Charge ; A History of the Lanai ey Aeronautical 
Laboratory . 1917-1 958 (in press, NASA, Washington, D.C.), plus 

projected outline of its sequel: 


* * * 


Foreword (by Dr. Donald Hearth, director of the Office of Space 
Science and Technology, University of Colorado, and former director 
of the NASA Langley Research Center) 

I n troduc 1 1 on 

1 . Foundat i ons 

2. Langley Personality, Formative Years 

3. The Oar i abl e-Dens i ty Wind Tunnel 

4. With a View to Practical Solutions 

5. The Cowling Story: Experimental Impasse and Beyond 

6. The Challenge of Teamwork 

7. The Priorities of World War II 

8. Exploring Unknown Technology: The Case of Jet Propulsion 

9. The Transonic Problem 

10. Defining the Research Airplane 

11. The Slotted Throat and Area Rule 

12. Hypersonics and the Transition to Space 

Ep i 1 ogue 

Notes 

Appendi xes 
Bibliographic Essay 


Projected outline of sequel: 


13. 

A New Mission 


14. 

Genesis of the 

Lunar Orbit Rendezvous Concept 

15. 

What Happened 

to the First "A" in NASA 0 

16. 

Planning for an American SST 

17. 

Managing Lunar 

Orbi ter 

18. 

The Hypersonic 

Research Engine Project 

19. 

A Viking Raid 


PI us, 

updated appendixes 

, etc . 
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COMPUTATION OF TRANSONIC VORTEX FLOWS 
PAST DELTA WINGS-INTEGRAL EQUATION APPROACH 


Osama A. Kanchl 

Professor of Mechanical Engineering 
and Mechanics 
Old Dominion University 
Norfolk, Virginia 23508 

Abstract 

The steady ful 1 -potential equation is written in the form of Poisson's 
equation, and the solution of the velocity field is expressed in terms of 
an integral equation. The solution consists of a surface integral of 
vorticity distribution on the wing and its free-vortex sheets and a volume 
integral of source distribution within a volume around the wing and its free- 
vortex sheets. The solution is obtained through sucessive iteration cycles. 

The source distribution is computed by using a mixed finite-difference scheme 
of the Murman-Cole type. The method is applied to delta wings. Numerical 
examples show that a conical shock is captured on the suction side of the wing. 
It is attached to the lower surface of the leading-edge vortex but does not 
necessarily reach to the wing surface. 

Introduction 

Efficient aerodynamic design of fighter aircrafts is still a challenge for 
computational aerodynami cists due to the various flow regimes encountered during 
fighter aircraft missions such as during air-combat-fighter missions and air-to 
ground-strike missions. Throughout these missions, the flow field may vary from 
an essentially attached condition to a partly separated condition to a fully 
separated condition. In addition to these complex flow conditions, flow 
compressibility changes from low subsonic to supersonic Mach numbers during , 
take-off, cruise, combat action, manuevenng and landing. Bradley and Bhateley 1 
reviewed the status of computational methods with emphasis on fighter design 
applications. Their ratings of the existing computational techniques for 
fighter aircraft applications ranged from poor to fair, a status which urgently 
calls for more improvements and developments. 

In the present paper, we deal with the problem of transonic vortex flows 
around highly swept wings - a key aerodynamic problem for the future develop- 
ment of super maneuvering fighter aircrafts. The procedure is an extension of 
the method of reference 2 where we exploit the shock-capturing nature of the 
method through careful computation of the source distribution in the volume 
integral term which represents the full nonlinear compressibility in the flow- 
field. In addition, we discuss an alternative shock fitting technique which 
remains to be tried. 
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Concluding Remarks 


The numerical results of the integral-equation approach for the transonic 
vortex flow problem clearly show that the integral -equation solution is able 
to capture the shock. The captured shock is of curved shape, its foot is 
attached to the lower surface of the leading-edge-vortex sheet, does not 
necessarily reach to the leeward side of the wing and is located almost at the 

spanwise location of the leading-edge vortex core. In the cases considered 

here, the shock was found to disappear at almost ten percent of the root chord 
behind the trailing edge. The computed results of the compressible flow, as 
compared to those of the incompressible flows, show that the leading-edge 

vortex core moves a little inboard from the leading edge and a little downward 

nearer to the suction side. The present velocity-field formulation is more 
advantageous than the earlier velocity-potential formulation for both the shock 
and vortex-sheet calculations. 

It should be pointed out that no attempt has been made to show the convergence 
of the solution versus the number of grid points of the finite volume. More 
work is still needed to prove this point and moreover to sharpen the shock either 
by refining the grid in the shock region or by fitting the shock while using a 
relatively.course grid. In addition the MNDV-method should be replaced by the 
NHV-method 1 ’ , and an additional half turn of the leading-edge vortex system 
should be considered to cure the problem of having the foot of the shock attached 
to the feeding vortex cut in the forward cross-flow planes. Finally, the merits 
of monotone differencing should be considered relative to the classical type of 
differencing. In addition to these improvements, the code should be transferred 
to a parallel processor. 
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SIMPLIFIED FREE VIBRATION ANALYSIS OF 
COFS USING RECEPTANCE METHOD 


Rakes h K. Kapania 
Aerospace and Ocean Engineering 
Virginia Polytechnic Institute and State University 
Blacksburg, VA 24061 


Abstract 


Introduction : Over the last few years, NASA has been developing the space 

shuttle in order to place in orbit a new class of spacecraft called Large Space 
Structures (LSS). A flight research program, COFS**? (Control of Flexible 
Structures), has been proposed. This experiment will systematically evaluate 
algorithms and techniques for on-orbit response prediction, system identifica- 
tion and flexible body control of large space structures. The program involves 
a series of flights (see Fig. 1) building progressively from modeling and 
modal characteristics of LSS to the more complex issue of flexible-body interac- 
tive control. The fourth in the series of flights may involve a beam-like 
structure with a fairly large attached structure, such as an antenna as shown 
in Fig. 2. To increase the challenge of identification of closely space coupled 
modes, it is desirable to have a rapid way to determine stiffness and mass 
characteristics of the structural components that will yield a design with 
closely space modes. In this study, the receptance method-^, 4 (also known as 
mobility, or dynamic flexibility) is chosen to analyze the different design 
configurations. 

Receptance Method : The receptance method is based on an electrical analogy 

called "Force - Current" analogy that involves Kirchoff's first law, which 
states that in any network the algebraic sum of all the current flowing towards 
a node equals zero. Receptance is defined as the ratio of amplitude of the 
harmonic response of a point in a subsystem produced by a harmonic force to 
the amplitude of the force F. If R is the amplitude of the response due to 
F, then 


a(u>) = R (u>)/F (1) 

Receptance methods are conveniently used to solve complex problems. In them, the 
complex system is resolved into subsystems that can be analyzed rather easily. 
These subsystems are then recombined to find the response at a point or points 
in the original system. 

Analysis of COFS : The structure is divided into two subsystems, as shown in 

Fig. 3. Subsystem A consists of a cantilever beam having 6 degrees-of-free- 
dom at its free end (3 translations and 3 rotations). Similarly, subsystem B 
consists of a free-free beam and a rigid tip mass (to simulate antennas) 
at its one end. The tip mass is at a distance "t" from the beam end. This 
offset, that results in coupling axial and in-plane bending of the free-free 
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beam and also torsion and out-of-plane bending, is expected to play a major 
role in producing coupled modes. Both the subsystems are joined through 6 
degrees of freedom at the free end of the cantilever beam (see Fig. 4). The 
free vibration equation is given as 3 


[ a(u>) ] A + [ a(u>) ] B 


6x6 


= 0 


( 2 ) 


where [ a(u>) ] A and [ a(w) ] are the receptance matrices of the two subsystems. 

The natural frequencies of the system are obtained by solving Eqn. 2. A Bisection 
method is used to solve the above equation. It should be noted that the two 
receptance matrices contain various trigonometric and hyperbolic functions and 
all the frequencies can be obtained from the above 6x6 matrix. Finally 
efforts will be made to find small changes in w due to small modifications in 
system parameters. 
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El j = BENDIND RIGIDITY FOR IN PLANE VIBRATIONS 
EI 2 = BENDING RIGIDITY FOR OUT-OF-PLANE VIBRATIONS 

GJ = TORSIONAL RIGIDITY 

Lj = LENGTH OF VERTICAL BOOM 

L 2 = LENGTH OF THE ARM 

= OFFSET LENGTH 


Figure 3.- Two subsystems of the complex COFS experiment. 
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FREE VIBRATION ANALYSIS 

+l cL t fu >H - 0 


Receptance Matrix for Subsystem A 
r * Receptance Matrix for Subsystem B 


Figure 4.- Combination of two subsystems at the free end 
of the cantilever beam. 
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ELECTRON RADIATION EFFECTS ON A 
POLYARYLENE ETHER 


Richard L. Kiefer 
Department of Chemistry 
College of William and Mary 
Williamsburg, VA 23185 

Structures in space such as antennas and platforms will likely be constructed 
with polymer-matrix-composite materials because of their light-weight, high- 
strength, and low-thermal expansion. However in geosynchronous orbit, such a 
structure would receive a radiation dose of 10^ to 10™ rads over its 
necessary lifetime of 20 to 30 years. Since a dose of this level is known to 
initiate degradation and/or crosslinking in polymeric materials, a study of 
the effects of radiation on such materials has been under way in the Applied 
Materials Branch of the Materials Division. 

In this specific project, a novel polyarylene ether, synthesized by chemists 
in the Polymeric Materials Branch, was subjected to various doses of electron 
radiation. Chemical effects of the radiation were determined by infrared 
spectroscopy. The molecular structure of the polymer is 



The material, designated Polymer I in our study, is amorphous with a glass 
transition temperature (Tg) of 152°C. It is aTso saluble in several common 
solvents so that films could be easily made. Polymer I has some molecular 
features in common with polyether ether ketone (PEEK), a well -studied polymer 
used in composite materials. Unlike Polymer I, however, PEEK is largely 
crystalline and thus has limited solubility in organic solvents necessitating 
high-processing temperatures in any application. 

Thin films of Polymer I were cast on glass plates from a solution in 
N,N-dimethylacetamide (DMAc). It was found that very thin films (about 0.3 
mil) were necessary in order to obtain the best spectrum. Spectra were 
recorded from 4000 to 400 cm~^ with a Nicolet Model 3600 Fourier Transform 
Infrared Spectrometer (FT-IR) at a resolution of 2 cm“l. Spectra were 
recorded and stored before and after irradiation. Holes were cut in the films 
to match two pins on the sample holder. Thus, the films could be placed on 
the sample holder after irradiation in exactly the same position as before the 
irradiation to ensure that the same portion of the film was viewed by the 
infrared spectrometer. The difference between spectra taken before and after 
irradiation was used to determine the chemical effects. 
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Strips of film were irradiated in a vacuum chamber with 70 keV electrons. 
Irradiation times of 45, 88, 120, and 190 hours were employed to give a total 
dose ranging from 4.9 x 10^ to 2.1 x 10^0 rads. Two films were simultaneously 
irradiated each time. The dose rate was monitored with a Faraday cup and was 
constant throughout all irradiations at 1.096 x 10° rads per hour. 

The results showed a decrease in absorption for the irradiated films primarily 
at 1658, 1595, 1497, and 1244 cm"*. These peaks are associated with vibra- 
tions of the phenyl-ether and phenyl -carbonyl linkages. A decrease in absorp- 
tion after irradiation indicates that chain fission is occurring at both the 
ether and carbonyl bonds. A small decrease in absorption at 2968 cm -1 indicates 
that some methyl -hydrogen atoms are also removed during irradiation. These 
atoms apparently either form hydrogen gas or become attached to phenyl -carbon 
atoms but do not attach to oxygen atoms since there was no evidence for formation 
of hydroxyl groups. 

In the 190 hour irradiation, possible effects of oxygen on the irradiated film 
were tested. The films were transferred from the vacuum chamber in a nitrogen- 
filled bag to a nitrogen-filled glove bag. Both films were mounted on sample 
holders and transferred to the spectrometer under nitrogen. Spectra of both 
films were recorded and stored, then one film was exposed to air while the 
other remained under nitrogen. Several spectra of both films were recorded 
throughout the course of the day. There was no discernable difference between 
the spectra of the two films either as a function of time or as a function of 
exposure to oxygen. Thus, for Polymer I, the chemical effects shown by infrared 
spectroscopy are apparently occurring during the irradiation. 

The FT-IR data will be correlated with viscosity and gel permeation measure- 
ments performed on films irradiated simultaneously with the films used for 
infrared spectroscopy. In addition, mechanical measurements will be made on 
irradiated films of the same polymer. 
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DECONVOLUTION AND EQUI CORRELATION PROBLEMS 
OF LIGHTNING DATA PROCESSING 


Ali Kyra la 

Professor of Physics 
Arizona State University 
Tempe, Arizona 85287 

Abstract 

The problem of separating out the particular aircraft and collection system 
influence from the F106 Lightning mission data is approached by designing a 
digital filter to recover input data from the output data collected. For the 
first stage the input signal is to be constructed from the soliton solution for 
the lightning stroke derived in NASA Technical Memorandum 86349 (Jan. 1985). By 
using the observed data as output a causal filter is constructed. From this an 
inverse operator is constructed with due regard for causality and this operator 
is then used to determine input data from observed output data. The input data 
so determined should then be free of the influence of the particular aircraft 
and collection system. 

Since the filter is causa l >o nly past and present inputs and outputs are used 
to determine the system vector - ?. The model is strictly discrete with (N+l) 
inputs and outputs so that 



with the form of the matrix of input values given by 



94 



0 


0 o 


0 


-N 


..0 


The components of S are s_,s. 


,s N and those of Y are y.,y-i 


i uc uumpui ici 1 u i o aic a,,u *uz>c wi i a i c jr 9j — *i » • * • *y m* 

After use of the theoretical 1 solution as input data together with themeasured 
output data collected on actual flights an S matrix of the form 
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from which input data X= S * Y unaffected by the collection system can be 
determined. 

The second problem considered is the statistical problem of optimally 
combining multiple inputs of the same type of data (e.g. the time rates of change 
of electric displacement or magnetic induction) so as to arrive at a linear 
combination with minimal variance. For this the principle of equi correlation 
is used. 

It is shown that the optimal combination of the input data can be attained 
normalizing each channel against its own standard deviation and that the variance 
corresponding to this normalization is given in terms of the covariances by 
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TRAINING ON PROGRAMMABLE CONTROLLERS 


Ellis E. Lawrence 


Industrial Technology Department 
Elizabeth City State University 
Elizabeth City, North Carolina 27909 


ABSTRACT 

The programmable controller, PC, is a significant component in today’s 
industries and/or governmental agencies. It meets many of the industrial 
control needs today and has been designed to adapt to the changing needs of 
tomorrow. 

This report focuses on aspects of training, with special emphasis on 
programming in ladder logic, that should be considered before attempting to 
train individuals on PC's. A review of the literature, informal detailed 
discussions with several PC manufacturer personnel, and actually attending two 
PC manufacturer ( Allen-Bradley and Texas Instruments) training programs re- 
vealed that: (1) cost of training, (2) prerequisite for training, (3) use of 

simulators, (4) programming and (5) teaching aids (video based) are features 
that should be considered when developing a PC training program. A general 
instruction set m ladder logic was also developed in this study to assist PC 
users in their training. 
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LIMITS TO GROWTH: 


AN OPERATIONS STUDY OF THE SPACE STATION 


Lawrence Leemis 
Professor 

School of Industrial Engineering 
The University of Oklahoma 
Norman, Oklahoma 73019 


Abstract 

The projected traffic to and from the Space Station from ground thru orbital 
operations is evaluated by a simulation model. Included In the model are 
different Shuttle payload types (logistics, propellant, attached loads, etc.), 
a lower bound of 28 days for ground turnaround time, two launch pads at the 
Kennedy Space Center, propellant storage and bolloff losses at the station, 
one remote manipulating system at the station, and both Shuttle and orbital 
transfer vehicle manifesting procedures. Results are presented for both the 
transient and steady state cases with respect to system growth. 

Prior to the advent of reusable space vehicles, operational issues were of 
secondary importance to performance criteria. Operational aspects of the 
Space Station must be evaluated prior to design to insure an effective use of 
resources. The Space Station transportation system is divided into two 
subsystems: ground operations and Space Station operations. Ground 

operations involve Shuttle support, which includes taking payloads to and from 
the Space Station. Operations on the Space Station involve OTV (orbital 
transfer vehicle) and OMV (orbital maneuvering vehicle) missions. These two 
subsystems have propellant storage and RMS (remote manipulating system) usage 
in common. Propellant for the OTV and OMV is stored - in tanks at the Space 
Station. Shuttle missions are scheduled to deliver propellant to the Space 
Station based on the amount of propellant at the Station and the number of 
missions reguiring propellant. The RMS is used to unload payloads from the 
Shuttle and to load payloads onto the OTV and OMV. 

The model involves several classic problems In operations research. Queueing 
occurs as Shuttles wait for payloads to become available, and Space Station 
vehicles wait for payloads to arrive to the Space Station. Propellant 
delivery is an inventory problem involving probabilistic lead times and 
probabilistic demand. An optimization problem is also embedded within the 
system. Payloads waiting for integration into the cargo bay of any vehicle 
must be manifested subject to weight and space constraints imposed by the 
vehicle. Since the system is too complex for tractable analytic results, 
simulation is used to estimate measures of performance associated with the 
system. 

The model is implemented as a network/discrete event/continuous model using 
the SLAM simulation package. The network portion of the model involves 
approximately AO nodes and 40 nontrivial activities. Entitles within the 
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network represent Shuttles, payloads (logistics, propellant, attached 
payloads, low earth orbit missions, geosynchronous orbit missions and 
planetary missions), OTVs and OHVs. There are 15 discrete event subroutines 
in the model (approximately 500 lines of code). These subroutines perform 
functions such as manifesting payloads for Shuttle missions, boosting the 
station periodically so that it will remain in the same orbit, summarizing 
shuttle missions, scheduling propellant flights if the propellant inventory 
position is less than the reorder point, and manifesting payloads for OTV 
missions. The continuous portion of the model specifies differential 
equations which control propellant boiloff at the Space Station. The 
differential equations are solved numerically by the SLAM package. 

The output estimates performance measures associated with the system for 
various resource capacities. Included in the output are OMV utilization, OTV 
utilization, RMS utilization, Shuttle delay times, launch pad utilization, and 
payload delivery time. Spot charts (plots of the Shuttle or OMV flights over 
time) are plotted to indicate periods of congestion and inactivity. Payload 
demand and resource capacities are varied to assess their impact on the system. 
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RESOLUTION LIMITS FOR A HOLOCINEMATOGRAPHIC VELOCIMETER (HCV) 
FOR TURBULENT FLOW MEASUREMENTS 


J. A. Liburdy 

Associate Professor of Mechanical Engineering 
Clemson University 
Clemson, SC 29631 

ABSTRACT 

Turbulent flow is, in general, highly three-dimensional in nature, even 
in an exclusively two-dimensional mean flow. Predictive schemes require input 
as to the nature of the local turbulence in order to accurately model flow 
characteristics. Since there does not exist a closed form mathematical 
representation of turbulence, current computational procedures rely on 
empirical input to define their model. To this end,- a new method of 
experimental measurements is being developed which is to determine much of the 
turbulent structure of a wide range of fluid flows. 

This new technique, holocinematographic velocimetry (HCV), uses 

holographic imaging of a particle seeded flow. A high speed camera advances 
individual film frames which record a hologram of a portion of the flow. The 
hologram is then to be used to reconstruct images of the particles 
(approximately 40 gm diameter). Tracing particles from frame to frame will 
determine local velocity vectors. 

Advantages of holography over other recording techniques include its 
three-dimensional characteristics which can be used to fully determine a local 
velocity vector. However, the resolution of particle images in and out of a 
plane normal to the direction of viewing is not as good as within the normal 
plane. In order to obtain uniform resolution information two orthogonal views 
are to be obtained simultaneously. 

The design of a fluid flow experiment for HCV requires, among many other 
things, the proper seeding of the flow. The seed must (i) accurately follow 
the fluid fluctuations, (n) be readily discernable in the holographic image 
and (m) be of the proper number density. The first requirement necessitates 
small particles that are neutrally buoyant. The second requirement implies 
sufficiently large particles. The third requirement implies a sufficiently 
low density so as not to influence the flow or cause particle interaction, yet 
a high enough density to allow a detailed mapping to be obtained. 

The present study addresses some of the resolution limitations of a 
specific proposed holographic recording and reconstruction technique. An in- 
line Fraunhofer optical arrangement is to be used to record the holographic 
image. The recording distance is to be approxi amtely 80 far-field distances 
(as limited by the experimental configuration). Roughly, the expected 

resolution is to be ± 5 wm in the image plane. In the proposed HCV the 
reconstructed holograms will be scanned with a digital television camera and 
image analyzer system which is to locate particle centers. It is the lack of 
resolution in the image plane-normal direction which may significantly 
increase the data reduction effort. A given particle is expected to stay in 
approximate focus over 1-3 mm. Over this focal range the data are to be cross 
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referenced with the orthogonal view to ultimately obtain i 5 pm resolution in 
all three coordinates. Results of this study document the variation of image 
intensity versus deviations from the exact focal distance. An analytical 
study is also presented which determines the particle-to-background intensity 
ratio at the particle edge versus out of focus proximity of near-by particles. 

In order to experimentally determine the resolution limits in-line 
Fraunhofer holography was used to image 38 um size paricles which were vapor 
deposited onto a glass substrate. Particles were positioned relative to each 
other, both in and out of the focal plane. A series of reconstructed images 
were viewed through a television monitoring system with approximately a 3.0 mm 
by 3.0 mm field of view. Photographs were taken of the images on the 
television screen and studied with respect to the particle image intensity 
distribution. The individual particles were found to have a high intensity 
over approximately a ±4-6 mm span normal to the focal plane. This implies 
that rather restrictive limits will have to be placed on the threshold level 
of particle detection. By doing so the resolution normal to the focal plane 
can be enhanced to facilitate cross referencing of particles from orthogonal 
views. 


The analytical results were obtained by assuming a linear film response 
and superimposing the amplitude functions of two particles with a given field 
of view. The intensity distribution of the reconstructed image was then 
determined from the amplitude function. The edge contrast was calculated 
based on the definition of visibility. The results show that as the particles 
are separated normal to the focal plane there is a gradual contrast variation 
to about 25 particle diameters. 

One further resolution limit of the HCV system is the ability to 
reproduce the recording medium (film) configuration during reconstruction. 
Since rapid frame advancement is to be used in the hologram formation there is 
the. potential for film bending out of a plane normal to the viewing 
direction. This will result in locally different effective focal distances in 
the reconstruction process. An analytical study was carried out to determine 
the allowable magnitude of bending. This is based on the Rayleigh resolution 
criteria. For the HCV system under consideration the film would be allowed to 
bow approximately 1.2 mm out of plane and be within the Rayleigh resolution 
(approximately 0.5 ym). 
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A STABLE FACTORIZATION APPROACH TO LARGE SPACE STRUCTURE CONTROL 


John H. Lilly 


University of Kentucky 
Department of Electrical Engineering 
Lexington, KY 40506 


Abstract 


Over the past few years, a significant body of new theory has appeared 
dealing with the so-called stable factorization approach to design of 
compensators for linear MIMO systems (see, for example, [1] and the 
references contained therein). The key feature of the approach is that the 
ring of proper, stable transfer function matrices is in fact a principal 
ideal domain. Thus, it is always possible to factor any rational transfer 
matrix as a product ND* 1 (also rf^N) where N and D (N and D) are both 
proper, stable transfer matrices, and furthermore have a greatest common 
right (left) divisor of the identity matrix. This factorization of the 
plant P as P=ND~* (P=D~*N) Is known as a rlght-coprime factorization 
(r.c.f.) left-corprime factorization (l.c.f.) of P. 


Once an r.c.f. and l.c.f. of the plant are found, there are many results 
which become available to the designer in designing a compensation scheme 
for the plant. For example, once an r.c.f. and l.c.f. (and the 

corresponding Bezout factors) are found for P, one may: 


(1) Obtain a parametrization of all compensators which stabilize P. 

(2) Determine the compensator which minimizes the sensitivity of P to 
random disturbances, while at the same time stabilizing P. 

(3) Determine the compensator which maximizes the robustness to plant 
uncertainties, while at the same time stabilizing P. 

(4) Determine the compensator which simultaneously accomplishes (2) and 
(3). 


Thus far, the stable factorization theory has been applied to plants with 
unspecified structure (i.e. the plant transfer matrix is simply given as 
P) . However, if the plant structure is specified, more specific results may 
be obtained about the compensators in (l)-(4) above. 

Large space structures (LSS) are hyperbolic distributed parameter systems, 
and are commonly modeled for design purposes as lumped systems of high 
order. Furthermore, the LSS model always has a specific and well-known 
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fora. Thus, if an r.c.f and l.c.f. for this model can be obtained, then the 
stable factorization approach may be applied to the control of LSS, and the 
designer may determine compensators which have properties such as (l)-(4) 
above. 

The main result of this research thus far has been the determination of an 
r.c.f. and l.c.f. for an LSS together with the corresponding Bezout factors 
(using results in [2]). In fact, the factorization obtained is doubly 
coprime , and is in terms of known physical quantities of the LSS, i.e. 
rigid-body inertias, vibrational damping, and natural frequencies. Two 
immediate results for LSS are (1) a parametrization of all stabilizing 
compensators in terms of known physical quantities of the LSS, and (2) the 
fact that LSS are strongly stabilizable , i.e. can be stabilized using a 
stable compensator. 

With this factorization, it Is now possible to Investigate such important 
questions as sensitivity minimization and robustness enhancement by 
compensators which can be specified in terms of known physical quantities of 
the LSS. 
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CONTROL OF STRUCTURES WITH JOINT- TYPE NONLINEARITIES 
USING STOCHASTIC OPTIMIZATION 
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ABSTRACT 

This work is divided into two disjoint parts. Part I concerns the 

theoretical analysis of the structure. Part II considers the design of an 
actual damping device. 

Part I : 

Many space applications involve an antenna (or other large object) 
which requires a high degree of pointing accuracy. This pointing accuracy 
is achieved, in part, through the position control of the antenna’s 
supporting structures. This work is directed toward developing a control 
positioning system for these supporting structures. 

Because the structures planned for space use are very large, it has 
been suggested that they be collapsible. As a result these structures will 
have many joints; they are called joint dominated . Joint dominated 
structures consist of flexible struts connected by rigid joints with 
nonlinear stiffness and damping characteristics. A straightforward 
analysis of the complex behavior of joint dominated structures is not 
tractable. Here we propose to approximate the complicated nonlinear model 
by a siraplier stochostic model. In this work we concentrate only on the 
modeling aspects (but with an eye toward control design). 

Consider the simple model of a joint dominated structure shown in 
Figure 1, where the joints have the nonlinear stiffness characteristics 
shown in Figure 2. This system can be modeled by 


x - f(x), x(0) =• x Q , t>0 


( 1 ) 


where x(t)ex*R^ n . Because we have modeled the nonlinearity as piecewise 
linear, (1) can be modeled as a set of linear differential equations as 
follows: 


S^: f ( x) = A^x if x(t)e 
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Here, the sets l°°k like 


n 


i 




• • • 


x n 


j 

n-1 ,n 


where ^,-l+i are regions in the (x t> x 1+1 ) plane. 

With each point x(t) on a trajectory, we associate its corresponding 
dynamic model If we sample the trajectory, then we can associate each 

trajectory with a sequence of models [M n ] where If 

x(t n ) e Now we attempt to approximate the nonlinear model by 

developing stochastic models of the sequences [ M^] . The transition 

probabilities are derived from an a priori distribution on the state space. 
This distribution, essentially a design parameter, is chosen by considering 
several simulations of the actual nonlinear system. The stochastic model is 
validated through Monte Carlo simulations. 

Part II ; 

Large spacecraft are expected to be highly flexible and require systems 
that suppress vibrations excited following orientation maneuvers, 
articulation of subassemblies, and other necessary spacecraft operations. 
This paper presents the analysis and design of a passive vibration damping 
device. The device, the variable coefficient viscous inertial damper 
(VCVID), is self-contained with no relative motion between it and the 

structure that is damps. This allows the device to be sealed making it 
ideal for use in spacecraft. The principle of operation of the damper is 

that motions of the external case react via a viscous working fluid and 
cause a force between the case and a heavy mass called the proof mass 
(figure 3). The valve allows the coefficient of damping to be varied during 
various spacecraft operations. 

Here we discuss the performance of this device on a flexible 
structure. The main tool in the analysis is the Dual Generalized Hessenberg 
Representation (Dual GHR). The Dual GHR is a state space representation 

which exposes the input/output structure as reflected into the internal 
dynamics. In particular, the (Dual) GHR can be used to identify nearly 
uncontrollable and/or unobservable subspaces. Near unobservability can be 
Interpreted as near pole-zero cancellation. 

In this paper the dynamics of a combined flexible structure and damper 
are represented by its Dual GHR. The NASA Langley Research Center Grid 
Apparatus is used as the structure. The effect of various parameters of the 
damper on the stability of the combined structure/damper is discussed. 
Based on a preliminary analysis using a high-order finite element modal 

model of the grid, the grid model is shown to be nearly unobservable. This 

leads to a reduced-order model of the grid for designing the damper. The 
reduced-order model was selected for closed-loop design based on a specified 
controller structure. However, closed-loop performance evaluations are made 
using the full-order modal model. 
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CONTRIBUTIONS TO IDENTIFICATION AND CONTROL OF STRUCTURES 
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Abstract 


The research effort divides Into three separate but related problems which are 
expected to result in separate publications. One is a new algorithm to implement 
the ERA method of system identification which requires less computer storage 
and can enhance the accuracy of the results. The second Is the development of 
a confidence criterion for identified parameters. This not only indicates how 
good the results are but also is expected to help one decide the proper order 
of the system model. The third research area is to extend a slewing control 
method, developed and experimentally tested by Dr. Jer-Nan Juang, to slewing of 
two flexible members attached one to the other. This work is expected to be 
tested experimentally by Dr. Juang. It has applications to robotics and to 
deployment of flexible structures. 


Summary 


pr. Jer-Nan Juang and others at NASA Langley Research Center have developed an 
algorithm called ERA, or Eigensystem Realization Algorithm, for identification 
of structures from impulse response or from free decay data. The approach is 
an innovative generalization of the Ho-Kalman approach to identification, with 
new criteria developed to help determine system order and recognize noise 
contributions. ERA has been demonstrated to be a very effective structural 
identification method in application to modelling the Galileo spacecraft from 
test data. 

The first of the three contributions this summer is the development of a new ERA 
style algorithm based on Grahm- Schmidt orthogonal izati on of data vectors rather 
than on singular value decomposition. The main objective was to generate an 
approach which would not require use of a main-frame computer, but there are a 
number of important potential advantages to the approach which might improve 
the numerical accuracy of the results. 

The new approach allows one to build up to the proper order of the modal by pro 
ducing an orthogonal 1 zed sequence of data vectors, followed by a simplified 
identification step. The earlier approach required over-estimating the true 
order and then reducing to obtain the desired model. It has also been demonstra 
ted that the larger the over-estimate, the larger the error in the final identi- 
fied result due to measurement noise, which Is one aspect which may produce an 
accuracy advantage to the new approach. 
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The new approach also has a valuable characteristic allowing one to choose 
which columns of data in the Hankel matrix are most important to the identifica- 
tion. A significant improvement has been shown in identification results 
when the columns of the Hankel matrix are allowed to skip various steps. 

In fact, this ability to skip steps is an important practical improvement of 
ERA over the original Ho-Kalman formulation. However, ERA has had only ad hoc 
methods to determine which steps to skip. The new ERA algorithm gives a way to 
pick which columns contain the most independent information about the system 
dynamics, and therefore supplies a criterion for choosing which columns to 
include in the Hankel matrix. 

The new algorithm is presently being programmed and tested to determine its true 
numerical characteristics and to evaluate its advantages and disadvantages over 
the previous programs. 

The second subject of research is the generation of a mathematical expression 
for the influence of white Gaussian noise in measurement data on the parameter 
identification using ERA. This gives, for the first time, a criterion expressing 
the confidence level associated with the identified system parameter. The result 
relies on linearization of the relationship between the parameters and the data 
to obtain an approximate covariance of the modal frequencies and damping. 

In a paper by Dr. J. Juang and Mr. R. Pappa presented at the Astrodynamics 
Specialist Conference in August 1985 (Veil, Colorado), Monte Carlo studies were 
presented giving the scatter in identified parameters resulting from adding noise 
to simulated data. It was observed that when the model dimension is larger than 
the actual system dimension, the scatter in the parameters increases. This suggests 
that the mathematical expression for the covariance of the parameters can be used 
to help determine the proper system order to choose in the identification. 

The third area of research is to extend the method of slewing control of a flex- 
ible appendage developed by Dr. Juang to the slewing of two flexible members, 
one attached to the end of another. The control was experimentally verified for 
one appendage, and shown to be very effective. It uses strain gauges to give 
feedback informatijn on appendage vibration, which is particularly simple and 
inexpensive hardware. The control design is based on linear-quadratic control 
with frequency shaping. 

Problems where one encounters slewing of flexible members attached to each other 
include deployment of spacecraft, shape control of hinged satellites, and 
robotics. In order to increase the operating speeds of robots, many people 
suggest the use of lightweight composite materials, which will require control- 
lers that can supply good damping characteristics to vibration modes. So the 
results will be applicable to future generations of robots, as well as to space- 
craft. 
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The multiple flexible body slewing control problem is fundamentally nonlinear. 
The controller design will use a mathematical model to cancel nonlinearities in 
order to be able to use linear contoller design. The same strain guage feedback 
measurements will be used on each link. Experiments will be performed once the 
contoller is designed, in order to study the practicality of the approach and 
to evaluate the dependence of performance on the accuracy of the model, both 
through the interpretation of the feedback measurements and throuqht the 
cancellation of nonlinearities. 

Dr. Juang and co-workers will continue to work on each of these three research 
problems after my return to Columbia University in order to complete the work. 
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ABSOLUTE NONLINEAR ACOUSTIC MEASUREMENTS 
WITH LITHIUM NIOBATE TRANSDUCERS 


Larry Mattix 
Associate Professor 

Department of Chemistry and Physics 
Norfolk State University 
Norfolk, Virginia 


Recently, there has been considerable interest in the use of ultrasonic 
techniques to characterize composite materials and their precursors. This 
study is part of an on-going investigation which seeks to correlate the 
mechanical properties of polymeric materials, like ultrasonic attenuation, 
velocity, and nonlinearity, to curing parameters like viscosity, free volume, 
and glass transition temperatures. We seek to develop a technique for using 
contact transducers to measure material nonlinearity through harmonic 
generation by ultrasonic pulses. 


When a high intensity ultrasonic wave passes through a nonlinear medium 
the waveform is distorted. This distortion arises because of a deviation from 
Hooke's law, which relates the pressure and volume (or stress and strain) m 
the medium, and is manifested by the generation of harmonics of the ultrasonic 
driving frequency (ref. 1). The measurement of harmonically generated elastic 
waves in solids is a well established technique for characterizing material 
nonlinearity. The absolute amplitudes of the fundamental and the harmonics 
can be used to calculate the higher order elastic constants which appear in 
the nonlinear wave equation for solids. The solution of the wave 
equation to second order yields the following relationship between the 
amplitude ot the second harmonic Ay and the amplitude of the fundamental Ai 
(ref. 1). 


A 2 = A 2 (k 2 /4)(K J /K 2 )x 


( 1 ) 


where k is the propagation constant given by 2it/A in which X is the wavelength 
of the fundamental, represents the second- and third-order elastic 

constants corresponding to the crystalline propagation direction, and x is the 
distance in the material through which the fundamental wave has passed. Thus, 
a measurement of amplitudes A^ and A 2 as function of x will provide a measure 
of the material nonlinearity. 


Absolute measurements of this type are traditionally made using 
capacitive detectors, where the sample itself is one plate of a parallel plate 
capacitor. To achieve the required sensitivity the capacitor gap must be 
quite small and accurately calibrated. This would require both the sample 
face and the electrode of the detector to be lapped flat to within an optical 
wavelength, a process which would preclude the use of the technique for 
routine material characterizations. The goal of our study is to determine the 
feasibility of using acoustically bonded lithium niobate transducers to make 
these measurements. 


The samples used were 2 inch diameter polycrystalline aluminum cylinders 
of lengths 1.0, 1.5, 2.0, 2.5, and 3.0 inches. The ends of the cylinders were 



carefully machined to be parallel and flat. A special apparatus was designed 
to maintain the alignment of the transducers along the cylinder axis and to 
insure a good reproducible bond between them and the surface. The 
experimental apparatus used in this investigation are shown m the block 
diagram in the figure 1. The rf output from a synthesizer was gated by a 
timing generator with the amplitude of the resulting tone burst set by a 
programmable attenuator under computer control. This high intensity frequency 
burst (the fundamental at 4 megahertz) was used to drive the 4 MHz transducer. 
The output of the 8 MHz transducer (the harmonic response of the sample) was 
filtered, amplified and digitized, then transferred to the minicomputer for 
off-line storage and processing. 

Preliminary results show agreement with equation (1) with very little 
correction necessary for bond thickness and diffraction effects. Further 
measurements are now m progress. 
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THE DESIGN OF AN AIRFOIL FOR LOW 
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Mark D. Maughmer 
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An investigation into the design of an airfoil for a high-altitude, long- 
endurance remotely piloted vehicle (RPV) has been undertaken. The mission 
specified for this aircraft involves flight at altitudes of up to 70,000 feet 
and endurances of up to 120 hours. As revealed from the detailed set of re- 
quirements developed, the design of an airfoil tailored to the given mission 
is complicated by the high cruise lift coefficients dictated and the wide 
range of Reynolds numbers encountered. Particular difficulties are associated 
with the lower Reynolds numbers for which the transition process is character- 
ized by the presence of laminar separation bubbles. Furthermore, the weight 
due to the large amount of fuel on-board initially pushes the need for a fair- 
ly high maximum lift coefficient to allow reasonable take-off runs. Also, while 
some degree of laminar flow is desired. for low drag at cruise, it is neces- 
sary that the maximum lift coefficient is achieved with fully turbulent flow 
so that take-off performance is not sensitive to surface contamination of the 
airfoil. Because the penalty for even a small weight increase is significant 
when carried over such long durations, the design effort undertaken is for 
a flapless airfoil. This makes the achievement of high lift coefficients and 
low drag difficult; however, the potential payoff in simplicity and weight 
savings justifies the consideration of a flapless configuration. In this regard, 
the overall performance of a single-element profile is compared to a flapped 
section in which the possibilities of lower drag are traded-off against in- 
creased weight. 
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OPTIMAL: A LANGUAGE FOR MANIPULATING 

PARSE TREES 


Larry J. Morel 1 

Assistant Professor of Computer Science 
College of William and Mary 
Williamsburg, Virginia 23185 

Software tools for the rapid development of reliable software are of vital 
importance to the software engineer. These tools are so critical that it is 
not uncommon to find meta-tools, i.e., tools used for the production of other 
tools. One of these meta-tools, the parser generator, automatically generates 
a program that will parse input specified by a Backus-Naur Form (BNF) grammar. 
The parser generator accepts the grammar and produces a set of parse tables 
that describes the actions necessary to recognize the grammatical structure of 
input. A table-driven parser may then be instantiated with these tables, 
yielding a parser for a specific language. The grammar can be augmented with 
semantic actions to be performed as the structure of the input is recognized. 
These semantic actions may build an intermediate data structure to capture 
some information from the input for later processing. A complete data structure 
that represents the structure of the entire input is called a parse tree. For 
such a representation, tools can be constructed that are independent of the 
grammar at hand. Thus, parser generators are particularly useful for environ- 
ments in which BNF-describable information is accessed through a set of tools. 

Despite the fact that many programs have inputs that are describable by a 
BNF grammar, the use of parser generators has been limited to the production of 
programming languages processors. The complex interface to existing parser 
generators dissuades all but the expert for foraging into their realm. A user 
must know parsing theory and the particulars of the parsing algorithm for the 
parser generator to be used. Augmentation of the BNF grammar with viable 
semantics requires full knowledge of the grammar and may require modifying the 
grammar. Not all BNF grammars are acceptable to parser generators and slight 
changes can caused a grammar to be rejected. Even worse, the flow of control 
among the semantic actions is difficult to design, debug, and modify since 
parsing does not supply the traditional one-in/one-out control structure so 
desirable in high-level programming languages. The programmer is not prevented 
from directly accessing variables supplied by the parser and is sometimes even 
encouraqed to do so. The resultinq semantic actions are then inextricably 
intertwined with the parser. Since the programmar must structure information 
collected during parsing, the resulting data structure is frequently tailored 
to the task at hand, rather than building a parse tree, for instance. Such 
tailoring is efficient in computer resources but burdensome on programmer 
development and maintenance time. The result is a system that is difficult to 
maintain and code that is impossible to reuse. 

A new programming language, OPTIMAL, has been designed to abstract the 
interface to parser generators. OPTIMAL allows the programmer to build and 
manipulate parse trees without knowledge of the underlying parser, its parsing 
technique, or the BNF grammar that drives the parsing. OPTIMAL therefore pro- 
vides the power of a parser generator without the concomitant pitfalls. OPTIMAL 



requires no expertise, no complicated control flow to understand, no grammars 
to modify, and no intertwining of the task at hand with the parser internals. 

The programmer can then concentrate on producing a structured solution that is 
easier to maintain and modify. 

To hide the details of the parser and its parsing technique, the language 
is an implied preprocessing language. OPTIMAL handles the low-level detail of 
reading the input and translating it into a convenient form, in this case, a 
parse tree. OPTIMAL provides operations for manipulating the resulting parse 
tree. Since a parse tree can be produced from any parser, OPTIMAL programs are 
therefore independent of any particular parser. 

To hide the grammatical details from the programmer, the language provides 
a specification-by-example mechanism called a unit . To manipulate, say. Boolean 
expressions, the user defines 'boolexp' to be 'a < b'. The OPTIMAL processor 
then parses the supplied example ('a < b'), thus deducing the grammatical unit 
intended by the user. The unit 'boolexp' can then be used in conjunction with 
the parse tree operations supplied by OPTIMAL in order to find Boolean expressions 
in the input, instantiate Boolean expressions, and perform transformations 
involving Boolean expressions. 

The main difficulty in designing an OPTIMAL language processor involves 
parsing of units. A new form of parsing, called unit parsing has been developed. 
It has been shown that for an input I acceptable to a parser, any substring S 
of I can be parsed in time linearly proportional to the time spent parsing S 
while parsing I. Unit parsing has been implemented for the MYSTRO parser 
generation system and is being incorporated into an OPTIMAL language processor. 
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NON-LTE EFFECTS IN THE MESOSPHERE 
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Abstract 


A study has been made to determine how to do radiative transfer 
calculations in non-equilibrium atmospheres. The stimulus for the project 
concerned data received from the Limb Infrared Monitor of the Stratosphere 
(LIMS). In inverting these data to determine atmospheric parameters, 
equilibrium conditions are assumed throughout the atmosphere. Some of the 
data from the upper atmosphere were received from regions which are very 
likely not in equilibrium. An attempt has been made to determine what effect 
this might have on the parameters determined by LIMS. 


Summary 

The Limb Infrared Monitor of the Stratosphere (LIMS) is a six-channel 
limb-scanning radiometer which flew on the NIMBUS 7 spacecraft and operated 
from October 1978 until May 1979. It was designed to measure radiance in six 
spectral channels in the infrared. These measurements were inverted to obtain 
temperature and pressure as functions of altitude, and concentrations of 
ozone, water vapor, nitrogen dioxide, and nitric acid as functions of altitude 
in the stratosphere and lower mesosphere. A description of the instrument and 
various results are given in Reference 1 and the other references cited 
therein. The LIMS viewing geometry is shown in Figure. 1. 
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The inversion procedure mentioned above makes use of the radiative 
transfer equation. In applying this equation to the LIMS retrievals, the 
assumption is made that "local thermodynamic equilibrium" (LTE) applies to all 
parts of the atmosphere. This means that for any air parcel containing 
molecules in excited states, collisions between molecules are frequent enough 
to transfer energy and bring the parcel to an equilibrium state before the 
molecules have a chance to radiate their energy away. This assumption makes 
the radiative transfer calculation relatively straightforward, and is 
generally valid in the lower atmosphere where collisions are frequent due to 
high density. 

As altitudes increase into the mesosphere (50-85 km), atmospheric density 
becomes low enough so that LTE may no longer be a valid assumption. For a 
tangent height of 60 km for example (see Fig. 1), much of the radiation along 
the ray path comes from altitudes higher than 60 km, hence, the need for 
radiative transfer calculations without assuming LTE. Very little has been 
done in this area, and much of my time has been spent investigating possible 
approaches to the problem. 

A visit was made to the Air Force Geophysics Laboratory at Hanscom Air 
Force Base, Massachusetts. A group in the Infrared Technology Division there 
has developed a computer code for doing non-equilibrium radiative transfer 
calculations (2). Required for input for the program are kinetic temperature 
and number density profiles as functions of altitude, and also the so-called 
"vibrational temperature" as a function of altitude. This is a parameter 
which describes relative populations of excited states of molecules of a given 
species. It is an indication of how far the species deviates from equilibrium 
conditions within an atmospheric sample. By specifying a particular molecule, 
a spectral band, and various atmospheric profiles, the program calculates 
radiance along a ray path (see Fig. 1) as a function of tangent height. The 
group at the Air Force Geophysics Laboratory has provided us with a magnetic 
tape containing the program and various input files so that we may do sample 
radiance calculations with or without the assumption of LTE. The goal is to 
determine what effect non-equilibrium conditions might have on the results of 
the LIMS inversion procedure. Preliminary results indicate that the effect 
might be significant at a tangent height of 70 km, which is roughly the upper 
limit of the region of the atmosphere observed by LIMS. If the effect proves 
to be appreciable at lower altitudes as well, then the inversion technique 
might have to be modified. 
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ABSTRACT 


This project was conducted under the sponsorship of the Crew/Vehicle 
Interface Research Branch at the NASA Langley Research Center. One of the 
research objectives of this branch is to provide advanced enabling 
technologies and design methodologies to such areas as display 
components/subsystems and pilot input/output interfaces. These research 
systems are usually implemented on the branch's Advanced Display Evaluation 
Cockpit (ADEC) simulator. The ADEC utilizes the VAX-1 1/780 as a host 
computer and the Adage 3000 as a color graphics raster scan generator. 

The purpose of this project was to design and develop an advanced 
interactive control/displ ay interface system for monitoring transport 
aircraft engine/subsystems. This interface system also had to interact with 
the engi ne/subsystems during caution and warning situations. 

An important design factor to be considered is the trend to use more 
dedicated controls/switches in the cockpit of transport aircraft. In order 
to reverse this trend, this project was designed around three 
multi-function, (mu 1 1 i- I egend) , programmable display push-button (PDP) 
devices. Each PDP features an interactive display consisting of 560 pixels 
in a 16 X 35 array and a solid state Hall effect switch. 

The first phase of this project was to define the system logic which 
interfaces the ADEC simulator subsystem graphics with the PDP devices 
(controlled by the logic and refresh control unit) and the VAX-11/780 host 
computer. 

A system being controlled by multi-function devices dictates that the system 
be user-friendly and yet be sophisticated enough to assist aircraft crews 
during high workload periods. 

This was accompl ished by using branching logic during normal 
engine/ subsystem operation and tailored logic during caution and warning 
situations. User-friendly philosophy was stressed throughout the design 
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since some caution and warning systems in transport aircrafts are considered 
themselves to be a hazard. 

The final phase of this project was devoted to writing a Fortran program for 
the VAX-11/780 which supports the design objectives. 
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EPR CHARACTERIZATION OF METAL-CONTAINING EPOXIES 


Maria Pacheco 
Science and Math Division 
College of the Virgin Islands 
St. Thomas, U.S.V.I. 00802 


High performance polymers are being considered for use as 
films, coatings, adhesives and composite matrix resins in aerospace 
applications . 

Among these polymers, several types of thermoset epoxies 
are currently being used due to their extreme versatility, excellent 
processability, chemical resistance, low density and low cost. 

Since pure polymers have not been able to meet the requirements 
needed for all aerospace applications, research efforts have been 
directed towards the modification of these materials, one of these 
modifications being the incorporation of metal ions into these potential 
aerospace resins. 

Dramatic changes can occur in the properties of these 
polymeric materials when small quantities of metallic ions are 
introduced into the system, such as increases in mechanical strength, 
electrical conductivity and moisture resistance, although the 
reasons for these changes are not fully understood. 

In order to obtain more information on the electronic 
structure and coordination of the metal ions in these polymer 
matrices, the magnetic properties of the systems have been studied 
using Electron Paramagnetic Resonance (EPR) Spectroscopy. Studies 
were carried out for MY-720 epoxy matrices containing Fe, Co, Ni 
and Cr metal ions in 1:50, 1:25 and 1:10 concentration ratios, as 
well as for an Epon-828-Cr 1:5 matrix. 

Preliminary results indicate the existence of both high 
and low spin states for the metal ions in the resin, depending on 
the metal/resin formulation. For the Cr ion, similar occupancy sites 
in the polymer matrix have beep found for both epoxy systems studied. 

For the Fe ion, changes in the EPR signal as a function of metal 
concentration indicate the possible existence of more than one 
occupancy site and/or metal state in the polymer matrix. Further 
characterization studies and semiquantitative results will be 
discussed in detail. 



EXPERIMENTAL DETERMINATION OF THE 
STIFFNESS MATRICES FOR SYMMETRIC LAMINATES 


Howard V. L. Patrick 
Assistant Professor 
Aeronautical Engineering Department 
Embry-Riddle Aeronautical University 
Daytona Beach, Florida 32014 

Abstract 

Graphite-epoxy composite materials have become of great interest in the 
aeronautics and space industry because of their high strength and stiffness to 
weight ratios. Acoustically induced material fatigue is of concern because of 
the pressure fields generated by aircraft and spacecraft propulsion systems on 
structures • 


Scientists would like to mathematically predict the displacment response 
of composite panels to acoustically induced random pressure fields. The panel 
strains can be evaluated once the displacement response has been ascertained. 
Various investigators (1-3) have developed theories which predict the 
deflection response of symmetric laminate plates. Solution of these equations 
requires the determination of the extensional and bending stiffness matrices 
which can be predicted theoretically (4). This laminate theory is based upon 
the static experimental material properties of a single unidirectional lamina. 

Concern exists about the validity of using static material properties in 
a dynamic analytic model. The modulus of elasticity (a measure of stiffness) 
can be determined indirectly by mechanically vibrating cantilever beams (5-8). 


Classical lamination theory allows prediction of laminate stiffness 
characteristics where 
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l. <Vk <2 k ' Vl> 
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ij k k “k-1 
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D ij - 1/3 j, <Vk <4 - 4 -,> 


A^j is called the extensional stiffness, B^j the coupling stiffness, and 
Djj the bending stiffness. Knowing the transformed reduced lamina stiffness 

Qlj (4) and z^, which is the distance from the plate middle surface to the 
k-th layer, the stiffness matrices can be calculated. 

The force -moment relationships for a laminate, where e°, ey,°and 

y£y are the normal and shear strains at the middle surface, < xf and 
the curvature, are 
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where Nx, Ny, and N xy are the normal and shear forces per unit width, 

M x and My are the bending moments per unit width in the x and y directions 
respectfully, and is the twisting moment per unit width. 


For the symmetric case we are considering, the laminas are geometrically 
symmetric around the middle surface as well as materially similar; the 
coupling matrix vanishes. Further simplifications can exist if the 
laminate is extenslonally balanced, Aj2 = A26 - 0 , or balanced in bending 
where Djfc = D26 = 0. 


When performing experiments, the loads are applied and the resulting 
deformations are measured. It is, therefore, advantageous to invert the 
force-moment equations because the deformations are the dependent variables. 
For a symmetric laminate, the coupling matrix B^j vanishes and the inversion 
is simply (4): 
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Verification of one set of 


-1 f 

We are using the notation what • 


stiffness implies verification of the other set. 
readily be calculated from the inverted form. 


The stiffness matrices can 


The extensional matrix coefficients can readily be determined by applying 
unaxial forces to the laminate in preferred directions. Let us first consider 

the case where Ny = N X y = 0, then from Equ. 3 wa have 

o*o* o * 

e x = Ajj N x> £y ■ A^ N x and = A^ N x * We need only apply known uniaxial 

stress to a specimen, measure the surface strain, and determine the 
coefficients • 


Similarly we can apply a uniaxial force Ny where N x - N X y = 0 and 
measure the surface strains and determine 


A . = e°/N , 
12 x y 


A 00 = e°/N and 

22 y y 


and A_, = y° /N . 
26 xy y 


The stiffness coefficient A^g is obtained by placing a laminate 
specimen in pure shear. For this case, N x = N y - 0 and from Equ. 3 we 

We can also verify the shear coupling coefficients 


have A,, = y° /N 

66 xy xy 

where A,. = e°/N and A.. = e°/N 

16 x xy 26 y xy 

during the tension tests. 


both of which can also be determined 


Examination of the equation for bending, Equ. 4, shows that we can 
determine the bending stiffness coefficients by applying pure bending and 
twisting moments to the test specimen. For example let My = M xy = 0 we 

could then determine D n = k /M , = < /M and D,, = < /M . Using similar 

x x x x x z y x x o xy x 

techniques for the bending moment My, and the twisting moment Mxy, we can 

determine the remainder of the bending stiffness coefficients. 


Most of the bending stiffness coefficients can be found by using a simply 
loaded cantilever beam and measuring the strain near the support (6). We can 
find the twisting stiffness coefficient Dgg by twisting a square plate 
sample by placing forces perpendicular to the plate at the corners. This is 
accomplished by exerting equal forces on the first and third (diagonal) 
corners while other equal forces press on the remaining corners in an opposing 
direction. Jones (4) states that Dgg - PL^/4w; where P is the applied 
load at the comers, L is the diagonal (corner-to-corner ) distance of the 
plate and w is the vertical displacement at the corners with respect to the 
center of the plate. 

A few tension tests have been performed on 8-ply T300/5208 graphite-epoxy 
angle ply laminates. This laminate consisted of nominally 1-mm (0.040-in.) 
thick material having a lay-up of (0°, 45°, -45°, 90° ) s , s designated 
symmetric laminate. The results of these tensile tests are shown in Table I 

where orientation refers to the outer fiber alignment in the direction of the 
applied tensile loads. 
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Based upon a lamina thickness of 0*172 mm (0.005 in.) laminate theory 

i ? it 

predicts = 18.86-m/N and - -5.5756-Gm/N for a total 

laminate thickness of 1.016 mm (0.040 in.). Inspection of Table I shows that 

i ? 

A 11 ' A 22 withln 1,2 * and a 12 = A 21 within 0.9 %, we can, therefore, 

readily conclude that experimentally the coefficients are equal as predicted 
by the laminate theory. 

The main concern is to determine if the laminate theory is valid. Review 

V f f 

of Table I shows that A^ ( A 22^ was w^h* 11 approximately 3 % and A^ was 

within approximately 4 % when comparing laminate theory with experimental 
results. Based upon these results, we can readily conclude that the laminate 
theory is accurate, at least for the extensional case considered. 

It is recommended that the tension and shear tests be continued to 
determine the remaining extensional stiffness coefficients. The bending and 
twisting moment fixtures should be fabricated and appropriate tests performed 
to determine the bending stiffness matrix coefficients. The results of these 
tests will determine the validity of the laminate theory. It is further 
recommended that the vibrating cantilever beam tests be performed to determine 
if dynamic stiffness is different than for static conditions. These dynamic 
tests can be used for determining material damping which is required for the 
plate dynamic deflection response theory. 



TABLE I 


8-Ply Experimental Extensional Stiffness Data 


Sample 

Comment 

Orien . 
(deg) 

* 

(mm) 

A ii 

<°i> 

CM 

-«r b a z 

s-/ 

-*; 2 

CGf> 

<°f> 

E 11 

(GPa) 

E 22 

(GPa) 

V 12 

(unitless) 

A-l-T 

0 

1.054 

17.96 

— 

5.389 

— 

52.18 

— 

0.300 

A-2-T 

0 

1.062 

18.34 

— 

5.629 

— 

51.37 

— 

0.304 

A-3-T 

0 

1.044 

18.43 

— 

5.455 

— 

51.99 

— 

0.298 

A-4-T 

0 

1.034 

18.63 

— 

5.507 

— 

51.93 

— 

0.296 

A-5-T 

0 

1.064 

18.03 

— 

5.282 

— 

52.13 

— 

0.304 

A-6-T 

0 

1.186 

18.49 

— 

5.575 

— 

45.63 

— 

0.305 

A-7-T 

0 

1.036 

18.62 

— 

5.656 

— 

51.82 

— 

0.304 

B-l-T 

90 

1.034 

— 

18.36 

— 

5.511 

— 

52.70 

0.300 

B-2-T 

90 

1.029 

— 

18.45 

— 

5.647 

— 

52.70 

0.306 

B-3-T 

90 

1.039 

— 

18.79 

— 

5.356 

— 

51.24 

0.293 

B-4-T 

90 

1.034 

— 

17.93 

— 

5.511 

— 

53.96 

0.307 

B-5-T 

90 

1.041 

— 

17.39 

— 

5.611 

— 

52.22 

0.323 

B-6-T 

90 

1.034 

— 

18.00 

— 

5.670 

— 

53.73 

0.315 

mean 

— 

1.053 

18.36 

18.15 

5.499 


51.01 

52.76 

0.304 

6* 

— 

0.041 

0.27 

0.49 

0.134 


2.39 

1.00 

0.008 

% Diff. 

— 

3.7 

2.7 

3.7 

4.5 . 

3.6 

8.7 

5.6 

2.0 


ft » ? 

Theoretical Values : = 18.86 Gm/n, A^ = A^ = 05.756 Gm/n and t = 1.016 mm, 

DOD/NASA design data: E n = E 22 = 55.9 GPa and V = 0.31 

* Standard deviation 
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CONTROL OF FLEXIBLE BEAMS THROUGH ROOT ACTUATION 


Harry H. Robertshaw 
Professor 

Mechanical Engineering Department 
Virginia Polytechnic Institute and State University 
Blactsburg, Virginia 24061 


The problem of controlling flexible structures m space is rapidly 
becoming more than an analytical exercise. We are being faced with the 
prospect of having to actually perform vibration and position control 
on light, lightly damped structures in space. The investigation of root 
control of beams in the laboratory (on the ground) is uncovering many 
idiosyncrasies of these systems that the analyses do not predict. The 
research performed during the 1985 summer is addressing some of these 
problems. This abstract will cover the following areas: modelling and 

control overview, the problem of motor friction, controlling the beams 
in the gravitational field, and introduction of flexibility control to the 
base of the beam. 

We have addressed ourselves to the control of flexible beams via 
root actuation. We are attempting to do this task both in and out of 
the gravitational field. Our approach to this task is first to model 
the beams, actuators, and sensors and then add an analog linear state 
variable feedback law to place the eigenvalues of the resulting 

closed-loop system. The models for the actuator (a dc servo-motor 

with PWM power amplifier) and sensors (potentiometer for root angle and 
beam strain measurements) are made with standard assumptions. The 
beam continuum was discretized using a three— term Ritz approximation. 
The resulting model is an eighth— order ODE. Previous work (1984 
NA5A/ASEE Fellowship) demonstrated our ability to add significant 
damping to the system through state feedback using a geared motor 
with a standard dc amplif ler . Our present work with a direct-drive 
motor has shown that the motor friction is now a significant factor. 
The friction we encountered was approximately three-f ourths of the 
control torque required m one simulated slewing task. 

We have attempted to eliminate the effects of this significant 

friction using a ’’tuneable 11 open— loop bang— bang function. We have 
added an analog circuit to the control law that is an approximation of 
the static friction torque characteristics of the motor. The analog 
circuit used is an operational amplifier with two diodes in the feedback 
loop with the motor speed measurement (voltage) as the input. These 
diodes are placed in opposite directions; one of them serves to limit 
the maximum value of the amplifier output; the other serves to limit 
the lower value of this output. These maximum and minimum values are 
equal to the breakdown voltages of the diodes; therefore, the upper 
and lower values can be set independently by proper choice of the 
diodes. We did not have this asymmetric coulomb friction and used the 
same type diodes for each limit. We found that this circuit tended to 
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oscillate at zero speeds and added a low— pass filter to reduce the 
amplitudes of these oscillations- This circuit eliminated a major 
portion of the friction and allowed us to use the linear control law 
more effectively. 

In order to understand the effect of the gravitational field we 
are attempting to control the beams with the gravitational field 
parallel to the plane of rotation. The resulting equations of motion 
are nan-linear m the root angle (sine and cosine non— linearities). Our 
first approach to this control task has been to linearize the equations 
o-f motion at different root angles and then choose state feedback 
coefficients to set the eigenvalues of this linearized system. We 
wrote a one-mode approximation for the system and used a symbolic 
manipulation code (MACSYMA) to solve for the feedback gains as a 
function of the angle. 

Another approach taken to this control task was to investigate 
the possibility of accomplishing it by controlling the stiffness of the 
support the beam has at the root. This work has only begun. The 
preliminary work has involved the simulation of a simple one— degree— of— 
freedom spring-damper system with a variable spring rate. We have 
approached this problem via controlling an analog simulation (on an EAI 
2000) of this system by changing the stiffness externally. Our results 
have been mixed: some of the situations have demonstrated improved 
damping and others have not. More "experimental” work will be conducted 
along with an analysis of the time-variable parameter ODE that 
represents this variable flexibility situation. 
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REPRESENTING LINEAR SYSTEMS IN 
GRASSMANN MANIFOLDS 


George Rublein 
Department of Mathematics 
College of William and Mary 
Williamsburg, VA 23185 

Abstract 

We shall discuss some aspects of geometric representations available for linear 
control systems. A technique used by many authors exploits Grassmann manifolds 
as a context for describing many important invariants of these systems. We hope, 
eventually, to use this method to describe phase and gain margins for multi-input 
and multi -output systems. Such a scheme presumably requires the use of a 
"natural" metric on the manifold so that, for example, one may reasonably estimate 
separation of a Nyquist locus form crucial singularities. 

Meanwhile, we need to explain the variety of ways the Grassmann formalism is 
used. These may be briefly listed: 

1. For (complex) p-planes in C, we may use an mxp matrix of spanning vectors 
as 


C-(i) 

X 


where G is a rational proper transfer function and I is an identity matrix. This 
method is used by Brockett and Byrnes 1 and allows a reasonably straightforward 
description of the Nyquist locus. 

2. For p-planes in C, we may use a homogeneous description of the plane as 
an nxn non-singular matrix. By passing to an appropriate equivalence class, we 
can get a single valued version of the plane. This method is used by Saeks and 
Murrays forthe homogeneous representation of coprime (stable) factorizations of 
transfer functions. 
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PROGRAMMING OF ROBOT ARMS BY GEOMETRIC DESCRIPTION 


J. C. Sanwal 

Department of Mathematics 
College of William and Mary 
Williamsburg, Virginia 23185 


Abstract: 


We give a method for programming cooperating manipulators, which is guided 
by a geometric description of the task to be carried out. For this we must 
have a suitable language and a method for describing the workplace and the 
objects in it in geometric terms. A task level command language and its 
implementation for concurrently driven multiple robot arms is described. The 
language is suitable for driving a cell in which manipulators, end effectors 
and sensors are controlled by their own dedicated processors and these 
processors can communicate with each other through a communication network. A 
mechanism for keeping track of the history of the commands already executed 
allows the command language for the manipulators to be event driven. A frame 
based world modeling system is utilized to describe objects in the work 
environment and any relationships that hold between these objects. This system 
provides a versatile tool for managing information about the world model. 
Default actions normally needed are invoked when the data base is updated or 
accessed. Most of the first level error recovery is also invoked by the 
database by utilizing the concepts of demons. The package can be utilized to 
generate task level commands in a problem solver or a planner. 
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1.0 INTRODUCTION 


In designing a controller for a manipulator the manufacturer uses diverse 
schemes to control the motion of the manipulator arm. These can be thought of 
as machine instructions for the manipulator. We give a standard instruction 
set which can be constructed from the manufacturers set. Thus different 
manipulators can be integrated in a system in a uniform manner. This set is 
adequate to implement task level conditional commands for a manipulator when we 
interface them with a data base. 


2.0 STANDARD SET 


For a manipulator we can implement the following basic commands (ref. 1) 
Orlando): 

(INITIALIZE ARMNUMBER) — (1) 

(MOVE ARMNUMBER XYZOAT-LIST) — (2) 

(MOVE-ALONG ARMNUMBER DIRECTION-DISTANCE-LIST) — (3) 

(SET-ANGLES ARMNUMBER ANGLE-LIST) — (4) 

(POSITION ARMNUMBER) — (5) 

(STATUS ARMNUMBER IDENTIFICATION-NUMBER) — (6) 

(RESUME ARMNUMBER) — (7) 

(SUSPEND ARMNUMBER) — (8) 

(CANCEL ARMNUMBER) — (9) 


3.0 HISTORY MANAGEMENT 

If the receiving unit is enabled it puts on its command queue commands of 
type (l)-(4) and carries them out sequentially. Whenever a command has been 
completed the central processing unit moves it to its history stack along with 
the identification number given to it by the issuing unit. Notice that the 
issuers event count is used even though the device may maintain its own time 
stamp. Commands of the type (5) -(9) are executed immediately even if it is 
necessary to suspend a command of type ( 1 ) -( 4) . For commands of type (5) -(9) 
the commanding unit waits for the response. 

In a work cell system, as shown in figure 1, the connecting lines indicate 
the communication path between the central processing units. 


4.0 DATA BASE FOR WORLD MODELLING 

To be able to implement high level commands for the manipulators we have 
implemented a data base to represent the working environment using frames. In 
the implementation all the default actions normally needed are invoked when the 
data base is updated or accessed. Frames are structures to make declarative 
statements about objects and relationships between objects. Thus two 
typical frames in this world model, which describes an arm and an object of 
the system , would be: 
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(arml (type 
(base 
(init (value 
(located (value 
(using (hand 


(name (puma))) 

(radius (150.0))) 

((1064.4 681.461 1092.4 0.0 0.0 0.0)))) 
((1064.4 681.461 0.0 0.0 0.0 0.0)))) 
(handl)))) 


(a (type (value (cube))) 
(height (value (35.0))) 

(base (radius (24.75))) 

(on (object (p ))) 

(located (value ((1480.0 895.0 


35.0 3.14 0.0 0.0))))) 


5.0 CONCURRENT PROGRAMMING LANGUAGE 

Using the above history management scheme we can now implement a language 
to program multiple arm cell. Thus for the command 

(COND-PUT ARM OBJECT LOCATION SET-C0ND TEST-COND) — (10) 
all the information about the OBJECT is stored in a frame and suitable error 
recovery functions are invoked when the MOVE and PICK subcommands of the 
COND-PUT are issued. The COND-PUT waits for TEST-COND to be true before it 
starts executing and when completed sets the SET-COND to true. An interpreter 
for the language is implemented using the communication network, the time stamp 
and the history tables described above. 

To illustrate some of the commands and how they can be used to describe a 
task in the language we give a program to build a seesaw of figure 2. 

(PICK ARM1 BAR) ;These two are unconditional commands and thus 

(PLACE ARM1 BAR 0N-A) ;do not test or set conditions 

;Commands beqinnino with a C are conditional commands like (10) 

(CPICK ARM2 A2 C21 T) ; — (21) 

;Because the TEST-COND is true CPICK is not blocked 

;C21 will be set to T when CPICK has picked up the object A2 

(CPICK ARM1 A1 Cll T) ; — (11) 

; Cll will be set to T when the CPICK command has been executed 

(CM0VE ARM2 A2 RIGHT -END C22 C21) 

;Tests C21 from the CPICK command (21) above before continuing 
;C22 will be set to T when CMOVE has completed 
(CM0VE ARM1 A1 LEFT-END C12 Cll) 

;This command waits for (11) above to complete 
(CPLACE ARM1 A1 LEFT-END C14 (AND C12 C22)) 

;(AND C12 C22) guarantees that both objects A1 and A2 are immediately 
;above the bar so that if one cube is placed before the other the bar 
; would not fall 

(CPLACE ARM1 A2 RIGHT-END C23 C12) 

(CRETURN ARM1 H0ME1 Cl 3 C14) 

(CRETURN ARM2 H0ME2 C24 C23) 
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CONVENTIONAL CONTROL AND DATA DISPLAY MODIFICATIONS 
FOR ADVANCED MLS APPLICATIONS 


Charles Scanlon 
Professor 

Department of Mathematics, Computer Science, and Physics 
Arkansas State University 
State University, Arkansas 72467 


NASA Langley Research Center (LaRC) in a joint effort with the FAA, has begun a 
study of advanced applications of the microwave landing system (MLS). This 
program has the objective of defining an envelope of usable MLS approach paths 
considering pi lot -vehicle performance, pi lot -passenger acceptance, and ATC 
integration factors. 

The MLS consists of three major components: (1) ground-based scanning beam 

transmitters, (2) ground-based precision distance measuring equipment (DME/P), 
and (3) airborne MLS equipment for azimuth, elevation, and range 
determination. These components together with area navigation (RNAV) 
computation will be used to fly complex MLS procedures which include straight 
line segments and curved paths. The wide area coverage of the MLS permits 
curved path approaches representing a relatively new pattern or landing 
configuration for pilots, and requires modification of existing control and 
data display. 

Various levels of airborne equipment sophistication will be explored in the 
MLS study. The lower end of cockpit sophistication will be represented by 
conventional electromechanical instrumentation. The upper end of cockpit 
sophistication will be represented by an advanced configuration with fly-by- 
wire type controls and electronic displays. The conventional configuration 
studies will be conducted with the LaRC Visual Motion Simulator (VMS) and the 
LaRC DC9-30 fixed base full systems simulator, representing a category C 
aircraft. The advanced configuration studies will utilize the LaRC Transport 
Systems Research Vehicle (TSRV) fixed-base simulator, which also represents a 
category C aircraft. Studies in the VMS have begun. 

The VMS is a si x-degree-of -freedom, motion-base simulator capable of 
presenting realistic acceleration and attitude cues to the pilot. A general 
purpose, scientific mainframe computer with a nonlinear, high fidelity digital 
representation of DC9-30 twin-jet commercial transport airplane provides 

inputs to drive the VMS motion base system. Audio cues for engine thrust and 
aerodynamic buffet will also be provided. The simulator has a generic cockpit 
with conventional flight controls and instrumentation. Flight controls 

include a column and control wheel, rudder pedals, and throttle, speed brake, 
and flap controls located on a center console. Flight instrumentation 
includes conventional flight and navigation instruments and engine 
instrumentation. 

Initial tests in the VMS indicate that high pilot workload in complex MLS 
procedures will be a major limiting factor in designing usable approaches,. 

Efforts to reduce pilot workload have included the design and implementation 
of an autothrottle and the design of a flight director control panel. Pilot 

comments have also led to the design of a test to measure possible reduction 
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of pilot workload by displaying airplane track angle instead of magnetic 
heading on the horizontal situation indicator (HSI). 

The automatic throttle system consists of pilot interface hardware, aircraft 
state measurements, functional algorithms, and thrust setting commands. Pilot 
interface hardware consists of an autothrottle engage/disengage lever located 
on the throttle quadrant, disengage button and go-around palm switches on the 
throttle levers, commanded airspeed knob and readout located on the indicated 
airspeed instrument, and a fast/slow airspeed bug on the airplane attitude 
indicator. The algorithm for throttle setting is a modification of one used 
on the NASA TSRV and is based upon several variables including calibrated 
airspeed, commanded airspeed, along track acceleration, bank angle, and engine 
pressure ratio. 

The flight director control panel (FDCP) will be used to select guidance 
options available to the pilot in the VMS and subsequently in the fixed base 
DC-9 simulator. The FDCP is a modification of the FD-109 panel and will 

accommodate the original functions plus the addition of MLS and RNAV provided 
by the flight guidance and control system. The FDCP provides pitch select, 
mode select, and altitude hold inputs for the pilot. Selectable modes include 
GA (go-around), HDG (heading), NAV/LOC (VOR type navigation), APPR (ILS and 

MLS straight in type approach), and RNAV (area navigation including curved 
path MLS approaches). 

A plan of test has been developed to obtain comparative measures of the 
dependence of pilot workload and performance on displaying airplane track 
angle or magnetic heading on the HSI during a precision approach. Four to six 
airline pilots will be used as test subjects to fly airborne tracking 

navigation tasks in the VMS using the DC9-30 simulation software. 
Quantitative pilot workload and performance measurements will be recorded 
during each run and will include oculometer, heart rate, airplane track, 
airplane state, and pilot input measurements that will be time correlated for 
later analysis. Subjective evaluations will be obtained through pilot 
comments during and after each run and by pilot ratings, questionnaires, and 
debriefing after each session. 

The advanced application study is expected to continue through an approximate 
three year period. Joint simulation with the FAA using their ATCSF in 

conjunction with a NASA simulator and a FAATC-LaRC tie line is planned to 
address the MLS-ATC integration issues. 
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Now that the National Transonic Facility has been given operational status, 
experimenters are routinely listing among their test requirements the capa- 
bility of visualizing surface and freestream flow fields.* Unfortunately, 
there are presently no operational techniques for visualizing cryogenic flows. 
One candidate technique, vapor screening, is being examined in the present 
series of experiments. 

The vapor screen flow visualization technique, which is described in reference 
2, is being used for the first time to accomplish flow visualization in a cry- 
ogenic environment. In the experiments that have been performed to date, this 
technique has been used to visualize the flow over a 65° half-span delta-wing 
model at several angles of attack. The delta-wing model was mounted on a 
turntable in the test section in order to accomplish variations in the angle 
of attack. Fog was generated in the test section by running the tunnel near 
the freestream saturation temperature. The density of the fog was varied (by 
varying freestream temperature) in an effort to produce an optimum fog den- 
sity. The fog was illuminated by an intense sheet of light from a 15-mW 
helium-neon laser. Personnel from the Instrument Research Division designed 
the illumination system shown in figure 1. Video and still cameras, used to 
record the visualized flow field, were mounted inside a pod which was attached 
to one side of the test section. The swept (with respect to the tunnel walls) 
vertical light sheet (light source mounted on top tunnel wall) facilitated 
visualization of the three-dimensional vortical flow structure on the leeward 
side of the model. 

The vortex on the leeward side of the delta-wing model was visualized at 
several different tunnel conditions. Mach number (M) was varied from 0.4 to 

0. 8. total pressure (p ) was varied from 1.3 to 5.0 atmospheres; total temper- 
ature (T ) was varied t from 83 to 101 K and angle of attack was varied from 
13° to c 30°. Vapor screen photographs of the subject flow field are pre- 
sented in figure 2 at several angles of attack with M = 0.6, p = 3.0 atm. 
and T = 92.0 K. These flow visualization photographs are presently being 
interpreted relative to phenomenological implications. 

REFERENCES 

1. Hunter, W. W. and Foughner, J. T. (Editors), Flow Visualization and 
Laser Velocimetry for Wind Tunnels , NASA CR-2243, 1982. 

2. McGregor, I., "The Vapor Screen Method of Flow Visualization," Journal 
of Fluid Mechanics , Vol. 11, Pt. 4, Dec. 1961, pp. 481-511. 


137 










Figure 2 , View through D Window ot the flow over a 65 Half-Span 
Delta Wing at Various Angles of Attack (M-0.6. p ( “3,0 atm. 
and T— 92. OK) 
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Abstract 


A great deal of work has been dedicated to the development 
of structural control theories. Although the developments are 
extensive, the designer often finds it difficult to apply many of 
these theories to "real structural control problems." Indeed, it 
is of primary concern to bridqe the qap between engineering 
design and the structural control theories. 

The focus of the research into s t rue t ural con t ro 1 theories is 
diverse. Many of the researchers are concerned with practical 
implementation problems and toward that end promote decentralized 
control (Refs. 1, 2). Others promote a centralized modal control 
approach and toward that end point out that a control theory 
should not destroy certain ch a ract e r i s t i cs which are natural to a 
structure (Refs. 3-5). Still others, in search for the global 
optimum, are concerned with distributed controls (Refs. 6, 7). 

Much attention is also given to describing the robustness of the 
control theories in the presense of modeling errors, particularly 
in view of the fact that it is difficult to characterize 
structural stiffness in mathematical models (Refs. 8-10). 

All of the these approaches support the objective to 
uniformly dampen the motion of a spacecraft. As it turns out , a 
uniform damping control is a robust, decen t r a 1 i zed , natural 
control with near globally op t imal performance (R ef . llT i Thus , 
a uniform dampinq control answers the concerns raised in the 
previously sited references. 

In this paper, uniform dampinq is examined and specifically 
uniform damping is compared with globally optimal controls. This 
paper also investigates the effects of discretization of the 
controls in space, i.e., the effects of using a limited number of 
control forces. Also the effects of discretization of the 
controls in time are investigated. Finally, this work 
investigates the combined effects of discretization of the 
controls in time and space. As numerical examples, simple beams 
as well as complex truss structures are considered. 
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The purpose of this report is to define the procedures necessary to implement 
effective configuration control of electronically generated designs. As this 
writer was assigned to the Model Engineering Section (MoES), the work and 
techniques were used as examples for this report. However, the control mea- 
sures discussed can be applied beyond the generation of model designs to any 
of the varied design work done at Langley Research Center or throughout 
industry. 

Little previous research could be found concerning control of electronically 
generated designs. Even though the capability has existed to integrate de- 
sign and fabrication with computers for some time, few organizations have yet 
adopted such a system on a large scale, as the MoES has. Therefore, the 
approach taken was to first study and document the methods used in generating 
hand-drawn designs. Those steps were then updated to accommodate designs that 
were developed on the computer and machine-plotted. Finally, the total elec- 
tronic creation ana use of engineering designs were investigated. 

In the past designs have been formalized on engineering drawings so that the 
researchers, engineers, and technicians could visualize the end product. 

Because of the time consuming process of producing an original drawing the pro- 
cedure of requiring approval signatures on the original was sufficient to 
assure control. However, a new plotted drawing considered to be an original 
can be generated from an electronic definition in a matter of minutes. New 
plotted "originals" are being generated for even the smallest change or 
modification to a part. Though approvals are still required in this process, 
it is much easier to have multiple drawings carrying the same drawing number, 
but with minor differences in the configuration. At the least this leads to 
confusion, and more likely to higher costs and poor quality control. With the 
potential to produce and use design data without ever producing even one 
physical drawing, the control problem is compounded. 

Because of these problems effective configuration control procedures must be 
implemented. This paper presents some recommendations on where those controls 
should be placed and how they can be used without imposing unnecessary delays 
on the process of producing hardware from engineering designs. This paper 
also briefly discusses the controls used in the past with hand-generated 
original drawings, as well as those currently being used for computer-generated 
plotted drawings. 
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It is established that third-order velocity-product moments are of paramount 
importance in the faithful simulation of atmospheric turbulence. It is also 
established that turbulence superimposed upon a sheared mean-flow does not 
have the same statistical characteristics that a turbulence superimposed 
upon a constant mean-flow does. 

INTRODUCTION 

Faithful simulation of atmospheric turbulence is a problem of great concern to 
the aviation industry in general. In the early years of this discipline 
commonly used methods for artificially generating turbulence time-histories 

employed a Gaussian model of turbulence-*-, but in recent years the trend has 

been toward more realism with a modified Gaussian model^"^, where fourth-order 
moments are not obtained directly from the second-order moments; none of these 
two fundamental models, however, produces odd-order moments in the simulated 
"turbulence," and the net result is a random time-history that pilots who 
"fly" through this simulated disturbance in state-of-the-art fl ight 
simulators often describe as too "continuous" and too "monotonous" when 
compared to real atmospheric turbulence. Pilots refer to these collective 
features in the simulated time-histories as "lacking the [all-important] 
element of surprise" of authentic hydrodynamic turbulence. 

In view of this discrepancy between (as currently) simulated and real 
turbulence an effort was initiated to establish the role of third-order 
moments, as well as that of wi nd-shear, in the dynamical characteristics of 
the correlation structure of atmospheric turbulence; it is hoped that the 
results of this endeavor will stimulate the search for non-Gaussian turbulence 
models which duly reflect the true nature of the phenomenon, and 
simultaneously elevate the art of turbulence simulation to that of turbulence 
emulation . Briefly, the role of third-order moments is essentially to provide 
a "driving force" for the turbulence two-point velocity correlations, while 
the effects of wind shear are required in the simulation of turbulence 
occurring near the ground — as would be encountered in the situation where an 
aircraft is taking off or landing. Turbulence superimposed upon a wind shear 
of any type does not at all have the same statistical structure that a 
turbulence superimposed upon a uniform wind does; this striking and often 
unrecognized characteristic is due to the continual interplay between 
turbulence and the underlying mean-flow, an interplay which manifests itself 
through the nonlinear convective term in the Navier-Stokes equation. 


142 



CONCLUSIONS 


Third-order velocity-product moments are of paramount importance in the 
correlation structure of turbulence and should accordingly be included in any 
faithful numerical simulation of the same. These moments ingress as a "driving 
force" into the differential equation for the two-point velocity-correlation 
of the turbulence, and as such convey upon it its characteristically unique 
features; indeed without them, the dynamical aspects of turbulence reduce to 
those of Gaussian diffusion, and unjustly purloin from it its natural "self- 
interaction" as well as its interaction with the underlying mean-flow. 
Turbulence superimposed upon a sheared mean-flow does not have the same 
structure as does a turbulence superimposed upon a constant mean-flow; this 
distinction is due to the natural interplay between (real) turbulence and the 
mean-flow. The presence of shear produces an alteration of the "viscous- 
friction" factor in the turbulence, to the point where this effect is 
operative even in the "high" Reynolds number situation; in unsheared flow this 
is typically not the case, and viscous forces are negligible whenever the 
turbulence Reynolds number is "high." 
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The objective of this research is to investigate the implementation of "a 
posteriori” error estimates for adaptive refinements in finite element 
computations. For many model problems these methods have been well developed, 
yet they are still not in common use in practice. Here we discuss one such 
method, which has good prospects for wide use. 


We will consider elliptic partial differential equations which can be put into 
the following variational formulation: 


(1) Given f e H^, find u e such that 
B (u,v) = (f,v) for all v e H 2 , 

where Hj and H 2 are Hilbert spaces and B is continuous, coercive and bilinear form. 

Now, let and V 2 be finite dimensional subspaces of and H 2 respectively. 
Then the finite element formulation corresponding to (1) is : 


(2) Find ujj e Vj such that 

B ( u h> v ) = (f> v ) f° r a li v e v 2* 

The objective of adaptive methods is to design techniques to minimize the error, 
u-u^, measured in some appropriate norm. It has been shown by I. Babuska and 
W. C. Reinbolt that if J=('F,T,V 2 > is an admissible family, then 


with 


(3) Dj^ 11 u-u h 11 R <_ D 2 n 
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Sup B (u-u h , ^v) 
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where D^ and D 2 are bounded constants. 


We have applied these estimates to two model problems. In these applications 
the quantities provide a heuristic for optimization. To compute the^i’s one 
must solve a local auxiliary problem. The constants D, and D 2 are independent of 
the size of the elements used, but they do depend on tneir shape. 
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EVALUATION OF SUWER RESEARCH PROGRAM 
(NASA-ASEE) 

BY 

FELLOWS 

Please complete and return to Gene Goglia by July 26, 1985, NASA MAIL 
STOP 122. 

1. Were you able to obtain sufficient information regarding the program? 
YES NO (Circle One) 

Comments : 


2. Was the contact with your anticipated research associate an early one 
to allow some preparation for the research assignment? 

YES NO (Circle One) 

Comments: 


3. Were you given a choice of research topics? 
YES NO (Circle One) 

Comments : 


4. Have you found your research problem: 

A. Challenging? YES NO (Circle One) 
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PAGE TWO 


B. Within your chosen field of interest? YES NO (Circle One) 
Conments : 


5. Do you plan to embark on research of similar nature upon your return 
to your home institution? 

YES NO (Circle One) 

Comments : 


6. First- Year Fellows: Do you plan to return to Langley next year, if 

invited? 

YES NO (Circle One) 

Comments: 


7. Did you find the lectures that you attended of interest? Any sug- 
gestions for the future. 

YES NO (Circle One) 

Comments: 
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PAGE THREE 


8. Your comments regarding the stipend you received would be appreciated. 
Comments : 


9. Did you experience difficulty in locating suitable housing? 
YES NO (Circle One) 

Comments: 


10. Was the administration of the program to your satisfaction? 
YES NO (Circle One) 

Comments: 


11. Additional comments regarding the program as a whole. You may use the 
back of this page also. 
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PAGE FOUR 


12. If you anticipate publishing a paper as a result of your research this 
summer, please provide the title of paper and journal in which you 
anticipate it being published. 


13. If you plan on submitting a proposal for continuance of your research 
effort this summer, please provide title of proposal and agency to 
which proposal will be submitted. 


Signature 
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NASA-ASEE 


SUMNER FACULTY RESEARCH PROGRAM 
QUESTIONNAIRE FOR RESEARCH ASSOCIATES 

Please complete and return to Gene Goglia by July 26, 1985, NASA MAIL STOP 122. 

CIRCLE ONE 

1. Would you say that your Fellow was adequately prepared for YES NO 

his/her research assignment? 

Comments: 


2. Would you comment on the diligence, interest, and eithusiasm with which your 
Fellow approached his/her research assignment. 


3. Would you be interested in serving as a research associate YES NO 

again? 

Comments : 


4. Would you be interested in having your Fellow (if eligible) return a second 
year? 

Comments : 
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PAGE TWO 


5. Any recommendations regarding improvement of the program will be 
appreciated. 


Signature 
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